10、《地理空间数据处理与动态地图创建》

《地理空间数据处理与动态地图创建》

1. 地理坐标转换

在地理空间数据处理中，经常需要将地理坐标转换为像素坐标。以下是相应的代码实现：

yDist = geoMatrix[5]
rtnX = geoMatrix[2]
rtnY = geoMatrix[4]
pixel = int((x - ulX) / xDist)
line = int((y - ulY) / yDist)
return (pixel, line)

接下来是具体的操作步骤：
1. 打开源图像：

src = gdal.Open("/qgis_data/rasters/satimage.tif")

2. 获取地理变换对象：

geoTrans = src.GetGeoTransform()

3. 执行转换：

world2Pixel(geoTrans, -89.59486002580364, 30.510227817850406)

4. 验证输出结果是否为 (1296, 1346) 。

这种转换在小区域内