我有一些使用rsgislib分割图像的测试代码。从我运行的rsgislib文档中:
from rsgislib.segmentation import segutils
inputImg = 'clipped_geotiff_image.tif'
outputClumps = 'clipped_jers1palsar_stack_clumps_elim_final.tif'
outputMeanImg = 'clipped_jers1palsar_stack_clumps_elim_final_mean.tif'
segutils.runShepherdSegmentation(inputImg, outputClumps, outputMeanImg, minPxls=100)
在尝试通过以下方式显式添加坐标系和其他空间信息后,outputMeanImg不包含坐标系和其他空间信息:
rsgislib.imageutils.copySpatialAndProjFromImage(outputMeanImg, inputImg)
或者
rsgislib.imageutils.assignProj(outputMeanImg, rsgislib.imageutils.getWKTProjFromImage(outputMeanImg), None)
有关详细信息,我在下面提供了gdalinfo的输出。有人能建议如何设置相应的系统和其他投影信息,或将子数据集BAND1/数据转换为常规光栅文件吗?你知道吗
======
%gdalinfo剪辑\u jers1palsar \u堆栈\u束\u elim \u最终_梅恩·凯亚你知道吗
Driver: HDF5/Hierarchical Data Format Release 5 Files: clip_jers1palsar_stack_clumps_elim_final_mean.kea Size is 512, 512 Coordinate System is `' Metadata: BAND1_DATA_BLOCK_SIZE=256d
BAND1_DATA_CLASS=IMAGE BAND1_DATA_IMAGE_VERSION=1.2
BAND1_NO_DATA_VAL_NO_DATA_DEFINED=
BAND1_OVERVIEWS_OVERVIEW1_BLOCK_SIZE=256d
BAND1_OVERVIEWS_OVERVIEW1_CLASS=IMAGE
BAND1_OVERVIEWS_OVERVIEW1_IMAGE_VERSION=1.2
BAND1_OVERVIEWS_OVERVIEW2_BLOCK_SIZE=187d
BAND1_OVERVIEWS_OVERVIEW2_CLASS=IMAGE
BAND1_OVERVIEWS_OVERVIEW2_IMAGE_VERSION=1.2
BAND1_OVERVIEWS_OVERVIEW3_BLOCK_SIZE=93d
BAND1_OVERVIEWS_OVERVIEW3_CLASS=IMAGE
BAND1_OVERVIEWS_OVERVIEW3_IMAGE_VERSION=1.2
BAND1_OVERVIEWS_OVERVIEW4_BLOCK_SIZE=46d
BAND1_OVERVIEWS_OVERVIEW4_CLASS=IMAGE
BAND1_OVERVIEWS_OVERVIEW4_IMAGE_VERSION=1.2 Subdatasets:
SUBDATASET_1_NAME=HDF5:"clip_jers1palsar_stack_clumps_elim_final_mean.kea"://BAND1/ATT/DATA/FLOAT SUBDATASET_1_DESC=[256x1] //BAND1/ATT/DATA/FLOAT (64-bit floating-point)
SUBDATASET_2_NAME=HDF5:"clip_jers1palsar_stack_clumps_elim_final_mean.kea"://BAND1/DATA SUBDATASET_2_DESC=[2500x1500] //BAND1/DATA (16-bit unsigned integer)
SUBDATASET_3_NAME=HDF5:"clip_jers1palsar_stack_clumps_elim_final_mean.kea"://BAND1/OVERVIEWS/OVERVIEW1 SUBDATASET_3_DESC=[625x375] //BAND1/OVERVIEWS/OVERVIEW1 (16-bit unsigned integer)
SUBDATASET_4_NAME=HDF5:"clip_jers1palsar_stack_clumps_elim_final_mean.kea"://BAND1/OVERVIEWS/OVERVIEW2 SUBDATASET_4_DESC=[312x187] //BAND1/OVERVIEWS/OVERVIEW2 (16-bit unsigned integer)
SUBDATASET_5_NAME=HDF5:"clip_jers1palsar_stack_clumps_elim_final_mean.kea"://BAND1/OVERVIEWS/OVERVIEW3 SUBDATASET_5_DESC=[156x93] //BAND1/OVERVIEWS/OVERVIEW3 (16-bit unsigned integer)
SUBDATASET_6_NAME=HDF5:"clip_jers1palsar_stack_clumps_elim_final_mean.kea"://BAND1/OVERVIEWS/OVERVIEW4 SUBDATASET_6_DESC=[78x46] //BAND1/OVERVIEWS/OVERVIEW4 (16-bit unsigned integer) Corner Coordinates: Upper Left ( 0.0, 0.0) Lower Left ( 0.0, 512.0) Upper Right ( 512.0, 0.0) Lower Right ( 512.0, 512.0) Center ( 256.0, 256.0)
我找到了部分解决方案,但在rsgislib中没有。你知道吗
使用GDAL,我可以打开一个子数据集:
上面的代码为我感兴趣的子数据集设置了投影和geotransform,但是请注意,我将输出打开为只读。如果我将其设置为update,它将失败,并且它实际上会更新选定的数据集,即使它是以只读模式打开的。我认为这是一个bug,但它似乎是通过创建一个新的。辅助.xml文件。你知道吗
也许有一个更简单/更好的方法来做到这一点,但这是目前的工作。你知道吗
相关问题 更多 >
编程相关推荐