30GB超えの GeoTIFF ファイルをどう処理しようか TIPS

[mapconcierge]

とりあえず 39.21GB の GeoTIFF 中身チェック

• Size: 144,808 x 141,467 pixel (およそ 15万x15万ピクセル＝225億ピクセル)
• Origin = (139.460526322999897,35.653509161690003)
• Pixel Size = (0.000000470615000,-0.000000384235000) [deg], 4cm 分解能
• 139.460526322999897,35.653509161690003 → 139.460526793614897, 35.653509545925003

taichi@DRONEBIRD Extreme SSD % /Applications/QGIS3.20.app/Contents/MacOS/bin/gdalinfo inagi202109278ebx.tif 
Driver: GTiff/GeoTIFF
Files: inagi202109278ebx.tif
Size is 144808, 141467
Coordinate System is:
BOUNDCRS[
    SOURCECRS[
        GEOGCRS["WGS 84",
            DATUM["World Geodetic System 1984",
                ELLIPSOID["WGS 84",6378137,298.257223563,
                    LENGTHUNIT["metre",1]]],
            PRIMEM["Greenwich",0,
                ANGLEUNIT["degree",0.0174532925199433]],
            CS[ellipsoidal,2],
                AXIS["geodetic latitude (Lat)",north,
                    ORDER[1],
                    ANGLEUNIT["degree",0.0174532925199433]],
                AXIS["geodetic longitude (Lon)",east,
                    ORDER[2],
                    ANGLEUNIT["degree",0.0174532925199433]],
            ID["EPSG",4326]]],
    TARGETCRS[
        GEOGCRS["WGS 84",
            DATUM["World Geodetic System 1984",
                ELLIPSOID["WGS 84",6378137,298.257223563,
                    LENGTHUNIT["metre",1]]],
            PRIMEM["Greenwich",0,
                ANGLEUNIT["degree",0.0174532925199433]],
            CS[ellipsoidal,2],
                AXIS["latitude",north,
                    ORDER[1],
                    ANGLEUNIT["degree",0.0174532925199433]],
                AXIS["longitude",east,
                    ORDER[2],
                    ANGLEUNIT["degree",0.0174532925199433]],
            ID["EPSG",4326]]],
    ABRIDGEDTRANSFORMATION["Transformation from WGS 84 to WGS84",
        METHOD["Position Vector transformation (geog2D domain)",
            ID["EPSG",9606]],
        PARAMETER["X-axis translation",0,
            ID["EPSG",8605]],
        PARAMETER["Y-axis translation",0,
            ID["EPSG",8606]],
        PARAMETER["Z-axis translation",0,
            ID["EPSG",8607]],
        PARAMETER["X-axis rotation",0,
            ID["EPSG",8608]],
        PARAMETER["Y-axis rotation",0,
            ID["EPSG",8609]],
        PARAMETER["Z-axis rotation",0,
            ID["EPSG",8610]],
        PARAMETER["Scale difference",1,
            ID["EPSG",8611]]]]
Data axis to CRS axis mapping: 2,1
Origin = (139.460526322999897,35.653509161690003)
Pixel Size = (0.000000470615000,-0.000000384235000)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  COMPRESSION=LZW
  INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  ( 139.4605263,  35.6535092) (139d27'37.89"E, 35d39'12.63"N)
Lower Left  ( 139.4605263,  35.5991526) (139d27'37.89"E, 35d35'56.95"N)
Upper Right ( 139.5286751,  35.6535092) (139d31'43.23"E, 35d39'12.63"N)
Lower Right ( 139.5286751,  35.5991526) (139d31'43.23"E, 35d35'56.95"N)
Center      ( 139.4946007,  35.6263309) (139d29'40.56"E, 35d37'34.79"N)
Band 1 Block=256x256 Type=Byte, ColorInterp=Red
  Overviews: 72404x70734, 36202x35367, 18101x17684, 9051x8842, 4526x4421, 2263x2211, 1132x1106, 566x553, 283x277, 142x139
  Mask Flags: PER_DATASET ALPHA 
  Overviews of mask band: 72404x70734, 36202x35367, 18101x17684, 9051x8842, 4526x4421, 2263x2211, 1132x1106, 566x553, 283x277, 142x139
  Unit Type: metre
Band 2 Block=256x256 Type=Byte, ColorInterp=Green
  Overviews: 72404x70734, 36202x35367, 18101x17684, 9051x8842, 4526x4421, 2263x2211, 1132x1106, 566x553, 283x277, 142x139
  Mask Flags: PER_DATASET ALPHA 
  Overviews of mask band: 72404x70734, 36202x35367, 18101x17684, 9051x8842, 4526x4421, 2263x2211, 1132x1106, 566x553, 283x277, 142x139
  Unit Type: metre
Band 3 Block=256x256 Type=Byte, ColorInterp=Blue
  Overviews: 72404x70734, 36202x35367, 18101x17684, 9051x8842, 4526x4421, 2263x2211, 1132x1106, 566x553, 283x277, 142x139
  Mask Flags: PER_DATASET ALPHA 
  Overviews of mask band: 72404x70734, 36202x35367, 18101x17684, 9051x8842, 4526x4421, 2263x2211, 1132x1106, 566x553, 283x277, 142x139
  Unit Type: metre
Band 4 Block=256x256 Type=Byte, ColorInterp=Alpha
  Overviews: 72404x70734, 36202x35367, 18101x17684, 9051x8842, 4526x4421, 2263x2211, 1132x1106, 566x553, 283x277, 142x139
  Unit Type: metre

[mapconcierge]

思考メモ

• できる限りオープンソースのものだけで解決したい → MapTiler Desktop Pro には頼らない
• Google Earth Engine だと最大 10GB GeoTIFF までしかインポートできない
• とりあえず gdal_translate -projwin [west] [north] [east] [south] [input file name] [output file name] で適切なサイズに切り出すか...
• 扱いやすいように 10GB 未満に切り出すとすると 7万x7万ピクセル＝49億ピクセル程度で切り出す。

[mapconcierge]

{
  "type": "FeatureCollection",
  "features": [
    {
      "type": "Feature",
      "properties": {},
      "geometry": {
        "type": "LineString",
        "coordinates": [
          [
            139.460526322999897,35.653509161690003
          ],
          [
            139.460526793614897, 35.653509545925003
          ]
        ]
      }
    }
  ]
}

{
  "type": "FeatureCollection",
  "features": [
    {
      "type": "Feature",
      "properties": {},
      "geometry": {
        "type": "LineString",
        "coordinates": [
          [
            139.460526322999897,35.653509161690003
          ],
          [
            139.460526793614897, 35.653509161690003
          ]
        ]
      }
    }
  ]
}