GeoTIFF được tạo bằng GDAL không có phạm vi tối thiểu / tối đa khi được tải trong Arc / QGIS

Tôi đang tạo một số trình quét điểm nổi trong GDAL bằng trình điều khiển GeoTIFF. Khi tôi tải các hình ảnh thu được vào QGIS hoặc Arc, hệ thống ký hiệu mặc định đặt min ở -3,40282e38 và tối đa là 3,40282e38, do đó, raster trông chính xác màu xám. Có cách nào để tôi có thể viết phạm vi thực tế vào GeoTIFF trực tiếp để khi được tải vào chương trình GIS, nó sẽ tự động chia tỷ lệ biểu đồ độc đáo không? Tôi đã thử tạo một biểu đồ mặc định như thế này:

rasterMin, rasterMax = raster.GetRasterBand(1).ComputeRasterMinMax()
raster.GetRasterBand(1).SetDefaultHistogram(rasterMin, rasterMax, 255)

Nhưng tôi gặp một lỗi trông giống như một vấn đề ràng buộc:

  File "/usr/lib/python2.6/dist-packages/osgeo/gdal.py", line 846, in SetDefaultHistogram
    return _gdal.Band_SetDefaultHistogram(self, *args)
TypeError: not a sequence

Bất kỳ đề xuất nào về những gì tôi có thể làm, hoặc nếu tôi đã phạm sai lầm?

Nếu bạn đang tạo dữ liệu, thì bạn có thể theo dõi tối thiểu và tối đa khi đang di chuyển. Nếu bạn không có quyền kiểm soát trực tiếp dữ liệu nguồn, hãy nói nếu bạn đang kết hợp dữ liệu hiện có, thì việc sử dụng ComputeRasterMinMax(0)là ổn.

Dù bằng cách nào, sau khi bạn đã nhận được các giá trị tối thiểu và tối đa của mình, bạn nên gọi SetStatistics()trên băng tần và đảm bảo bạn đặt tập dữ liệu của mình Nonekhi bạn hoàn thành để đảm bảo mọi thứ được ghi vào đĩa.

Ngoài ra, một điều gì đó ở phía sau tâm trí tôi đang nói với tôi rằng việc sử dụng dataset.GetRasterBand(1).<some band function>sẽ không hoạt động chính xác do cách GDAL Python xử lý các tham chiếu băng tần. Tôi không thể nhớ các chi tiết, nhưng có một cái gì đó trong danh sách gửi thư của GDALDev một thời gian trước đây. Vì vậy, tốt nhất là tạo một biến băng rõ ràng : band = dataset.GetRasterBand(1), sau đó gọi band.SetStatistics(min, max, mean, stddev).

Tôi không biết cách sử dụng GDAL, nhưng bạn có thể sử dụng python để gọi công cụ 'tính toán thống kê' trong Arc.