ベクターのラスター化について
樋田泰之
QGIS version: 3.16.4-Hannover
QGIS code revision: 654e76b3cb
Qt version: 5.11.2
GDAL version: 3.1.4
GEOS version: 3.8.1-CAPI-1.13.3
PROJ version: Rel. 6.3.2, May 1st, 2020
プロセシングアルゴリズム...
アルゴリズム 'ベクタのラスタ化（rasterize）' を開始しています...
Input parameters:
{ 'BURN' : 0, 'DATA_TYPE' : 5, 'EXTENT' : '-22205.036800000,-21894.295200000,-104990.708300000,-104790.945700000 [EPSG:6675]', 'EXTRA' : '', 'FIELD' : 'field_4', 'HEIGHT' : 3000, 'INIT' : None, 'INPUT' : 'delimitedtext://file:///C:/Users/2142/Desktop/%E7%AD%89%E9%AB%98%E7%B7%9A%E7%94%A8210310.csv?type=csv&useHeader=No&maxFields=10000&detectTypes=yes&decimalPoint=,&xField=field_3&yField=field_2&zField=field_4&crs=EPSG:6675&spatialIndex=yes&subsetIndex=no&watchFile=no', 'INVERT' : False, 'NODATA' : -9999, 'OPTIONS' : '', 'OUTPUT' : 'C:/Users/2142/Desktop/aaa.tif', 'UNITS' : 1, 'WIDTH' : 3000 }
GDAL command:
gdal_rasterize -l None -a field_4 -tr 3000.0 3000.0 -a_nodata -9999.0 -te -22205.0368 -104990.7083 -21894.2952 -104790.9457 -ot Float32 -of GTiff "file:///C:/Users/2142/Desktop/%E7%AD%89%E9%AB%98%E7%B7%9A%E7%94%A8210310.csv?type=csv&useHeader=No&maxFields=10000&detectTypes=yes&decimalPoint=,&xField=field_3&yField=field_2&zField=field_4&crs=EPSG:6675&spatialIndex=yes&subsetIndex=no&watchFile=no" C:/Users/2142/Desktop/aaa.tif
GDAL command output:
ERROR 4: `file:///C:/Users/2142/Desktop/%E7%AD%89%E9%AB%98%E7%B7%9A%E7%94%A8210310.csv?type=csv&useHeader=No&maxFields=10000&detectTypes=yes&decimalPoint=,&xField=field_3&yField=field_2&zField=field_4&crs=EPSG:6675&spatialIndex=yes&subsetIndex=no&watchFile=no' does not exist in the file system, and is not recognized as a supported dataset name.
Execution completed in 0.09 seconds
Results:
{'OUTPUT': 'C:/Users/2142/Desktop/aaa.tif'}
出力レイヤの読み込み
次のレイヤは正しく生成されませんでした。
• C:/Users/2142/Desktop/aaa.tif
QGISメインウィンドウの"ログメッセージパネル"をチェックすると、アルゴリズムの実行に関する詳細情報が表示されます。
--------------------------------------------------------------------------------------------------
メインのログは、
--------------------------------------------------------------------------------------------------
processing.run("gdal:rasterize", {'INPUT':'delimitedtext://file:///C:/Users/2142/Desktop/%E7%AD%89%E9%AB%98%E7%B7%9A%E7%94%A8210310.csv?type=csv&useHeader=No&maxFields=10000&detectTypes=yes&decimalPoint=,&xField=field_3&yField=field_2&zField=field_4&crs=EPSG:6675&spatialIndex=yes&subsetIndex=no&watchFile=no','FIELD':'field_4','BURN':0,'UNITS':1,'WIDTH':3000,'HEIGHT':3000,'EXTENT':'-22205.036800000,-21894.295200000,-104990.708300000,-104790.945700000 [EPSG:6675]','NODATA':-9999,'OPTIONS':'','DATA_TYPE':5,'INIT':None,'INVERT':False,'EXTRA':'','OUTPUT':'C:/Users/2142/Desktop/aaa.tif'})
--------------------------------------------------------------------------------------------------
恐らくデータベースに起因するのかと思い
ファイル名が日本語入っているためかと思って、半角英数にしても同様の症状です。
CSVテキストレイヤの追加の際に、Shift-JISを選択しても変わりませんでした。
他に一度SHPデータで保存してやってみても×
SHPデータにした後に、”属性のリファタリング”で「小数点付数値」を「64bit整数値｣にしても不可でした。
因みにSHPデータで一度保存してから、ラスター化しようとすると
------------------------------------------------------------------------------------------------------------
gdal_rasterize -l aaa -a field_4 -tr 3000.0 3000.0 -a_nodata 0.0 -te -22205.0368 -104990.7083 -21894.2952 -104790.9457 -ot Float32 -of GTiff aaa.shp C:/Users/2142/Desktop/aaaa.tif
GDAL command output:
ERROR 1: Attempt to create 0x0 dataset is illegal,sizes must be larger than zero.
ERROR 1: Cannot create C:/Users/2142/Desktop/aaaa.tif
Execution completed in 0.10 seconds
Results:
{'OUTPUT': 'C:/Users/2142/Desktop/aaaa.tif'}
出力レイヤの読み込み
次のレイヤは正しく生成されませんでした。
• C:/Users/2142/Desktop/aaaa.tif
QGISメインウィンドウの"ログメッセージパネル"をチェックすると、アルゴリズムの実行に関する詳細情報が表示されます。
----------------------------------------------------------------------------------------------------------------------------
と挙動が違うので、SHPからラスタ化したほうが解決策は近いのかもしれません。
何か解決策などを御教示頂けると助かります。
adachi
樋田泰之
adachiさま