Python中的Geotiff图像绘图

发布于 2025-02-09 19:22:27 字数 5033 浏览 0 评论 0原文

我是Python中的新的Geotif图像。我有两个同一区域的Geotiff图像。但是,在绘制这两个图像之后,它们不匹配。以下是两个图像的元信息。

图像1

Files: sli_hh_ML_I_221.tif
Size is 640, 7368
Coordinate System is:
PROJCRS["unknown",
    BASEGEOGCRS["unknown",
        DATUM["unknown",
            ELLIPSOID["unknown",1737400,0,
                LENGTHUNIT["metre",1,
                    ID["EPSG",9001]]]],
        PRIMEM["Reference meridian",0,
            ANGLEUNIT["degree",0.0174532925199433,
                ID["EPSG",9122]]]],
    CONVERSION["Polar Stereographic (variant A)",
        METHOD["Polar Stereographic (variant A)",
            ID["EPSG",9810]],
        PARAMETER["Latitude of natural origin",90,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8801]],
        PARAMETER["Longitude of natural origin",0,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8802]],
        PARAMETER["Scale factor at natural origin",1,
            SCALEUNIT["unity",1],
            ID["EPSG",8805]],
        PARAMETER["False easting",0,
            LENGTHUNIT["metre",1],
            ID["EPSG",8806]],
        PARAMETER["False northing",0,
            LENGTHUNIT["metre",1],
            ID["EPSG",8807]]],
    CS[Cartesian,2],
        AXIS["(E)",south,
            MERIDIAN[90,
                ANGLEUNIT["degree",0.0174532925199433,
                    ID["EPSG",9122]]],
            ORDER[1],
            LENGTHUNIT["metre",1]],
        AXIS["(N)",south,
            MERIDIAN[180,
                ANGLEUNIT["degree",0.0174532925199433,
                    ID["EPSG",9122]]],
            ORDER[2],
            LENGTHUNIT["metre",1]]]
Data axis to CRS axis mapping: 1,2
Origin = (-56469.605829227773938,5384.527309082705870)
Pixel Size = (89.627161160045119,-21.283702822737492)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  -56469.606,    5384.527) ( 95d26'48.61"W, 88d 7'46.10"N)
Lower Left  (  -56469.606, -151433.795) ( 20d27' 1.68"W, 84d40'26.23"N)
Upper Right (     891.777,    5384.527) (170d35'46.03"E, 89d49'12.04"N)
Lower Right (     891.777, -151433.795) (  0d20'14.66"E, 85d 0'32.78"N)
Center      (  -27788.914,  -73024.634) ( 20d50' 2.06"W, 87d25'25.54"N)
Band 1 Block=640x3 Type=Float32, ColorInterp=Gray
  NoData Value=nan

图像2

*ch2_sar_ncxl_20191017t234531221_d_sli_xx_fp_hh_d18_wrp.tif
Driver: GTiff/GeoTIFF
Files: ch2_sar_ncxl_20191017t234531221_d_sli_xx_fp_hh_d18_wrp.tif
Size is 2294, 6273
Coordinate System is:
PROJCRS["unknown",
    BASEGEOGCRS["unknown",
        DATUM["unknown",
            ELLIPSOID["unknown",1737400,0,
                LENGTHUNIT["metre",1,
                    ID["EPSG",9001]]]],
        PRIMEM["Reference meridian",0,
            ANGLEUNIT["degree",0.0174532925199433,
                ID["EPSG",9122]]]],
    CONVERSION["Polar Stereographic (variant A)",
        METHOD["Polar Stereographic (variant A)",
            ID["EPSG",9810]],
        PARAMETER["Latitude of natural origin",90,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8801]],
        PARAMETER["Longitude of natural origin",0,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8802]],
        PARAMETER["Scale factor at natural origin",1,
            SCALEUNIT["unity",1],
            ID["EPSG",8805]],
        PARAMETER["False easting",0,
            LENGTHUNIT["metre",1],
            ID["EPSG",8806]],
        PARAMETER["False northing",0,
            LENGTHUNIT["metre",1],
            ID["EPSG",8807]]],
    CS[Cartesian,2],
        AXIS["(E)",south,
            MERIDIAN[90,
                ANGLEUNIT["degree",0.0174532925199433,
                    ID["EPSG",9122]]],
            ORDER[1],
            LENGTHUNIT["metre",1]],
        AXIS["(N)",south,
            MERIDIAN[180,
                ANGLEUNIT["degree",0.0174532925199433,
                    ID["EPSG",9122]]],
            ORDER[2],
            LENGTHUNIT["metre",1]]]
Data axis to CRS axis mapping: 1,2
Origin = (-56469.605829227773938,5384.527309082705870)
Pixel Size = (25.000000000000000,-25.000000000000000)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (  -56469.606,    5384.527) ( 95d26'48.61"W, 88d 7'46.10"N)
Lower Left  (  -56469.606, -151440.473) ( 20d26'58.70"W, 84d40'25.49"N)
Upper Right (     880.394,    5384.527) (170d42'50.59"E, 89d49'12.26"N)
Lower Right (     880.394, -151440.473) (  0d19'59.10"E, 85d 0'32.00"N)
Center      (  -27794.606,  -73027.973) ( 20d50'12.96"W, 87d25'24.93"N)
Band 1 Block=2294x1 Type=Float32, ColorInterp=Gray
  NoData Value=nan
Band 2 Block=2294x1 Type=Float32, ColorInterp=Undefined
  NoData Value=nan*

I am new geotif image ploting in python. I have two Geotiff images of the same region. However, after plotting these two images they are not matching. The below is the meta information for both the images.

Image 1

enter image description here

Files: sli_hh_ML_I_221.tif
Size is 640, 7368
Coordinate System is:
PROJCRS["unknown",
    BASEGEOGCRS["unknown",
        DATUM["unknown",
            ELLIPSOID["unknown",1737400,0,
                LENGTHUNIT["metre",1,
                    ID["EPSG",9001]]]],
        PRIMEM["Reference meridian",0,
            ANGLEUNIT["degree",0.0174532925199433,
                ID["EPSG",9122]]]],
    CONVERSION["Polar Stereographic (variant A)",
        METHOD["Polar Stereographic (variant A)",
            ID["EPSG",9810]],
        PARAMETER["Latitude of natural origin",90,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8801]],
        PARAMETER["Longitude of natural origin",0,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8802]],
        PARAMETER["Scale factor at natural origin",1,
            SCALEUNIT["unity",1],
            ID["EPSG",8805]],
        PARAMETER["False easting",0,
            LENGTHUNIT["metre",1],
            ID["EPSG",8806]],
        PARAMETER["False northing",0,
            LENGTHUNIT["metre",1],
            ID["EPSG",8807]]],
    CS[Cartesian,2],
        AXIS["(E)",south,
            MERIDIAN[90,
                ANGLEUNIT["degree",0.0174532925199433,
                    ID["EPSG",9122]]],
            ORDER[1],
            LENGTHUNIT["metre",1]],
        AXIS["(N)",south,
            MERIDIAN[180,
                ANGLEUNIT["degree",0.0174532925199433,
                    ID["EPSG",9122]]],
            ORDER[2],
            LENGTHUNIT["metre",1]]]
Data axis to CRS axis mapping: 1,2
Origin = (-56469.605829227773938,5384.527309082705870)
Pixel Size = (89.627161160045119,-21.283702822737492)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  -56469.606,    5384.527) ( 95d26'48.61"W, 88d 7'46.10"N)
Lower Left  (  -56469.606, -151433.795) ( 20d27' 1.68"W, 84d40'26.23"N)
Upper Right (     891.777,    5384.527) (170d35'46.03"E, 89d49'12.04"N)
Lower Right (     891.777, -151433.795) (  0d20'14.66"E, 85d 0'32.78"N)
Center      (  -27788.914,  -73024.634) ( 20d50' 2.06"W, 87d25'25.54"N)
Band 1 Block=640x3 Type=Float32, ColorInterp=Gray
  NoData Value=nan

Image 2

enter image description here

*ch2_sar_ncxl_20191017t234531221_d_sli_xx_fp_hh_d18_wrp.tif
Driver: GTiff/GeoTIFF
Files: ch2_sar_ncxl_20191017t234531221_d_sli_xx_fp_hh_d18_wrp.tif
Size is 2294, 6273
Coordinate System is:
PROJCRS["unknown",
    BASEGEOGCRS["unknown",
        DATUM["unknown",
            ELLIPSOID["unknown",1737400,0,
                LENGTHUNIT["metre",1,
                    ID["EPSG",9001]]]],
        PRIMEM["Reference meridian",0,
            ANGLEUNIT["degree",0.0174532925199433,
                ID["EPSG",9122]]]],
    CONVERSION["Polar Stereographic (variant A)",
        METHOD["Polar Stereographic (variant A)",
            ID["EPSG",9810]],
        PARAMETER["Latitude of natural origin",90,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8801]],
        PARAMETER["Longitude of natural origin",0,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8802]],
        PARAMETER["Scale factor at natural origin",1,
            SCALEUNIT["unity",1],
            ID["EPSG",8805]],
        PARAMETER["False easting",0,
            LENGTHUNIT["metre",1],
            ID["EPSG",8806]],
        PARAMETER["False northing",0,
            LENGTHUNIT["metre",1],
            ID["EPSG",8807]]],
    CS[Cartesian,2],
        AXIS["(E)",south,
            MERIDIAN[90,
                ANGLEUNIT["degree",0.0174532925199433,
                    ID["EPSG",9122]]],
            ORDER[1],
            LENGTHUNIT["metre",1]],
        AXIS["(N)",south,
            MERIDIAN[180,
                ANGLEUNIT["degree",0.0174532925199433,
                    ID["EPSG",9122]]],
            ORDER[2],
            LENGTHUNIT["metre",1]]]
Data axis to CRS axis mapping: 1,2
Origin = (-56469.605829227773938,5384.527309082705870)
Pixel Size = (25.000000000000000,-25.000000000000000)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (  -56469.606,    5384.527) ( 95d26'48.61"W, 88d 7'46.10"N)
Lower Left  (  -56469.606, -151440.473) ( 20d26'58.70"W, 84d40'25.49"N)
Upper Right (     880.394,    5384.527) (170d42'50.59"E, 89d49'12.26"N)
Lower Right (     880.394, -151440.473) (  0d19'59.10"E, 85d 0'32.00"N)
Center      (  -27794.606,  -73027.973) ( 20d50'12.96"W, 87d25'24.93"N)
Band 1 Block=2294x1 Type=Float32, ColorInterp=Gray
  NoData Value=nan
Band 2 Block=2294x1 Type=Float32, ColorInterp=Undefined
  NoData Value=nan*

I tried with several methods but am unable to find why these two images do not display same.

如果你对这篇内容有疑问,欢迎到本站社区发帖提问 参与讨论,获取更多帮助,或者扫码二维码加入 Web 技术交流群。

扫码二维码加入Web技术交流群

发布评论

需要 登录 才能够评论, 你可以免费 注册 一个本站的账号。
列表为空,暂无数据
我们使用 Cookies 和其他技术来定制您的体验包括您的登录状态等。通过阅读我们的 隐私政策 了解更多相关信息。 单击 接受 或继续使用网站,即表示您同意使用 Cookies 和您的相关数据。
原文