PostGIS导入导出栅格数据

最新推荐文章于 2025-05-17 23:33:16 发布
原创 最新推荐文章于 2025-05-17 23:33:16 发布 · 1.7w 阅读 · 34 ·
CC 4.0 BY-SA版权
本文为博主原创文章，转载请注明原文出处！
文章标签：

#postgresql #PostGIS #GIS #数据库

本文详细介绍了如何在PostGIS中导入和导出空间栅格数据。使用raster2pgsql工具导入全球影像数据，通过psql查看数据，利用PostGIS内置函数如ST_AsGDALRaster、ST_AsPNG等导出数据。由于PostGIS默认不开启GDAL支持，需要修改配置文件启用。最后，通过Python的Psycopg库连接PostgreSQL并导出TIFF图像文件，虽然导出的图像丢失了colormap，但像素值保持正确。

摘要生成于 C知道 ，由 DeepSeek-R1 满血版支持， 前往体验 >

上一篇博文PostGIS导入导出ESRI Shapefile数据介绍了如何导入空间矢量数据到PostgreSQL中，紧接上一篇，本文将介绍如何使用PostGIS导入导出空间栅格数据。
实验环境和上一篇的博文一样： Ubuntu 16.04 LTS 。
数据为全球影像（ESRI ArcGIS提供的示例数据），下载连接：全球影像百度网盘下载
可以使用GDAL的gdalinfo命令查看其详细信息：
gdalinfo
在QGIS中查看如下：
wsiearth.tif

PostGIS提供了raster2pgsql工具用于栅格数据的导入（可以使用man raster2pgsql命令查看帮助文档）。
使用如下命令进行影像数据的插入，具体参数的含义这里不在累赘，上篇文章有说明。关于raster2pgsql的参数可以使用man命令进行查看。


                

                
                
        

    

    
  


        

            

                      
                        最低0.47元/天 解锁文章
                        

                          
200万优质内容无限畅学

                          
                        

                      
            

        


    



                



	

		

			

				
					
postgis 栅格数据_PostGIS导入导出栅格数据

				
			

			

				

					weixin_39606361的博客
				

				

					12-19
					
					725
					
				

			

		

		

			
				
上一篇博文PostGIS导入导出ESRI Shapefile数据介绍了如何导入空间矢量数据到PostgreSQL中，紧接上一篇，本文将介绍如何使用PostGIS导入导出空间栅格数据。实验环境和上一篇的博文一样： Ubuntu 16.04 LTS 。数据为全球影像(ESRI ArcGIS提供的示例数据)，下载连接：全球影像百度网盘下载可以使用GDAL的gdalinfo命令查看其详细信息： 在QGIS...

			
		

	



                

            
              



	

		

			

				
					
PostGIS批量导入栅格数据

				
			

			

				

					给永远比拿愉快
				

				

					02-14
					
					9441
					
				

			

		

		

			
				
如果我们有一批以文件存储的影像数据如何利用PostGIS批量的导出到PostgreSQL数据库中进行管理呢？ 
（单个数据的导入参见我上篇博文：PostGIS导入导出栅格数据）我的实验环境如下：  
OS: Ubuntu 16.04 LTS  
PostgreSQL:9.5.5 
PostGIS: 2.2 我的实验数据是全球气象数据（不过，我只下载了中国部分19区的平均气温数据（TIF格式）做示例演

			
		

	



              


  
              

                


	

		

			

				
					
PostGIS实现栅格数据入库【raster2pgsql】

					
最新发布

				
			

			

				

					Gerry-wu
				

				

					05-17
					
					1050
					
				

			

		

		

			
				
raster2pgsql是PostGIS提供的命令行工具，用于将GDAL支持的栅格格式（如GeoTIFF、JPEG、PNG等）导入PostgreSQL数据库，支持批量加载、分块切片、创建空间索引及金字塔概览，是栅格数据入库的核心工具。raster2pgsql通过灵活的参数配置和高效的批量处理能力，成为PostGIS栅格数据导入的首选工具。在实际使用中，需根据数据规模、格式特性及业务需求，合理选择分块策略、索引方式和事务模式，以平衡导入速度与数据完整性。更多开源GIS相关知识技能分享学习交流。

			
		

	





	

		

			

				
					
PostGIS 栅格导入导出

				
			

			

				

					苦逼码农的博客
				

				

					08-11
					
					1231
					
				

			

		

		

			
				
【代码】PostGIS 栅格导入导出。

			
		

	



		

			
		



	

		

			

				
					
PostGIS导入栅格切片

				
			

			

				

					qq_52231465的博客
				

				

					11-17
					
					699
					
				

			

		

		

			
				
其中，-s指定柳格教据的投影信息，表示创建空间索引，-C表示据，-M表示使用多人数据库连接，F表示通制转换为首通栅格，指定栅格据的块大小，public.raster table是要创建的表名。其中，-f指定要执行的SQL文件路径。

			
		

	





	

		

			

				
					
Postgis使用工具raster2pgsql导入栅格数据

				
			

			

				

					万里归来少年心
				

				

					03-14
					
					4423
					
				

			

		

		

			
				
1.栅格数据

从网站下载天气数据。



2.使用工具gdainfo查询栅格文件信息


>gdalinfo tmax1.bil
Driver: EHdr/ESRI .hdr Labelled
Files: tmax1.bil
       tmax1.hdr
Size is 2160, 900
Coordinate System is:
GEOGCRS["WGS 84",
    DATUM["World Geodetic System 1984",
        ELLIPS...

			
		

	





	

		

			

				
					
PostGIS导入导出ESRIShapefile数据

				
			

			

				

					02-23
				

			

		

		

			
				
这里介绍如何导入我们常用的ESRIShapefile数据到PostgreSQL数据库中，我们可以使用PostGIS提供的shp2pgsql和pgsql2shp工具进行导入和导出操作，还可以使用GDAL库提供的ogr2ogr工具，ogr2ogr工具支持更加多样的数据...

			
		

	





	

		

			

				
					
Postgis使用工具raster2pgsql批量导入栅格数据（一）

				
			

			

				

					万里归来少年心
				

				

					03-15
					
					4260
					
				

			

		

		

			
				
1.栅格数据

WorldClim是一个高空间分辨率的全球天气和气候数据的数据库。可以从该网站下载天气数据。网址如下：https://www.worldclim.org。

下面是12个月份温度统计数据。



2.使用工具raster2pgsql 将栅格数据导入到postgis

首先，生成sql文件。


...

			
		

	





	

		

			

				
					
workshop-postgis-raster:讲解并举例说明PostGIS栅格用法的研讨会

				
			

			

				

					05-11
				

			

		

		

			
				
研讨会PostGIS栅格
该研讨会旨在解释和举例说明PostGIS栅格的用法。
创建和还原数据库
为了开始研讨会，我们将加载一个现有的数据库。
 请使用pgAdmin创建一个新的空数据库，然后将现有的转储postgis_raster.backup从该存储库还原到新创建的新数据库中。
数据库结构
恢复的数据库由以下结构组成：
 schema_name （您应该使用您的名称重命名）
 上市
栅格
向量
铁路（线）
 地点（点）
 porto_parishes（多边形）
 Schema schema_name和栅格为空，您将在进行练习时向其中添加表。
 表格vectors.porto_parishes具有葡萄牙波尔图地区存在的所有教区。
 请先浏览数据库，然后再继续。
加载栅格
我们将从将栅格文件Landsat8_L1TP_RGBN.tif和srtm_1arc_v3.tif加载到栅格模式开始。 

			
		

	





	

		

			

				
					
基于PostGIS的高分一号卫星影像管理

				
			

			

				

					09-22
				

			

		

		

			
				
基于PostGIS的高分一号卫星影像管理与Web应用研究,专业论文，供参考

			
		

	





	

		

			

				
					
postgresql+postgis安装和空间数据的导入

				
			

			

				

					09-07
				

			

		

		

			
				
postgresql+postgis安装和空间数据的导入

			
		

	





	

		

			

				
					
Postgis使用工具raster2pgsql批量导入栅格数据（二）

				
			

			

				

					万里归来少年心
				

				

					03-15
					
					2565
					
				

			

		

		

			
				
另外一种批量导入栅格数据的方式，是使用工具gdalbuildvrt生成vrt文件 。

1.栅格数据

WorldClim是一个高空间分辨率的全球天气和气候数据的数据库。可以从该网站下载天气数据。网址如下：https://www.worldclim.org。

下面是12个月份温度统计数据



2.生成rvt文件


gdalbuildvrt -separate tmax_multi.vrt tmax*.bil

验证vrt文件的正确性


gdalinfo tmax_multi.vrt...

			
		

	





	

		

			

				
					
PostGIS导入矢量数据和栅格数据

				
			

			

				

					不去想结果，一直在路上
				

				

					12-19
					
					5529
					
				

			

		

		

			
				
PostGIS如何批量导入矢量数据和栅格数据1.矢量数据1.1 shp格式1.2 其它格式2.栅格数据
1.矢量数据
1.1 shp格式
PostGIS提供了导入shp数据的工具shp2pgsql，参数如下：




参数
描述




-s<srid>
设置数据的坐标，比如4326


-d|a|c|p
-d 删除原表并重新创建;-a 追加shp数据到当前表后;-c 创建一个新表，默...

			
		

	





	

		

			

				
					
postgresql导入栅格数据

				
			

			

				

					plusplusfei的博客
				

				

					03-03
					
					4368
					
				

			

		

		

			
				
1、首先，需要为数据库加载postgis的相关拓展。（https://postgis.net/install/）
– Enable PostGIS (as of 3.0 contains just geometry/geography)
CREATE EXTENSION postgis;
– enable raster support (for 3+)
CREATE EXTENSION postgis_raster;
– Enable Topology
CREATE EXTENSION postgis_to

			
		

	





	

		

			

				
					
gdal读取postgis数据库中的栅格数据（python）

				
			

			

				

					fusion_cell的博客
				

				

					07-29
					
					2084
					
				

			

		

		

			
				
文章目录前言一、向postgis导入栅格二、gdal读取postgis栅格

前言
gdal读取tif格式的栅格影像非常方便，它还支持读取postgis的栅格影像。不过，一开始我以为读取方式与tif会有较大差异，所以就上网搜索gdal如何读取postgis栅格。很不幸，内网几乎没有，外网没查过，一是不会魔法上网，二是硬看英文太累。功夫不负有心人，最后还是找到了利用gdal读取postgis栅格的JAVA实现（https://blog.csdn.net/qq_24309981/article/detail..

			
		

	





	

		

			

				
					
postsql导入tiff文件

				
			

			

				

					m0_58510078的博客
				

				

					03-04
					
					1021
					
				

			

		

		

			
				
满屏的乱码，执行完之后，数据库中也没有产生对应的表，那一般出现乱码的常见原因是编码问题。可是tiff文件没有编码，经过查找资料推断原因大概为解译tiff二进制时出了问题，建议可以先生成一个sql文件，然后在执行sql文件，将数据插入。只启用postgis是不够的，处理栅格数据，还要启用postgis_raster，之后重新运行sql文件，或者最开始使用的指令，都可以将tiff导入数据库。因为目前不用这个参数，数据也没啥问题，所以并没有研究为啥使用-k参数会出现这个错误，有知道的伙伴可以评论区告诉我。

			
		

	





	

		

			

				
					
栅格数据怎么做成csv表

				
			

			

				

					03-31
				

			

		

		

			
				
<think>好的，我现在需要帮助用户解决如何将栅格数据转换为CSV文件的问题。首先，用户提到了使用ArcGIS和PostGIS的方法，我需要先回顾一下这些工具的具体步骤。

首先，用户提供的引用中，引用[2]提到了ArcGIS中面转栅格，再栅格转点，最后导出CSV。这可能涉及到几个步骤：创建面状shp文件，转换为栅格，提取栅格中的点坐标，然后导出为CSV。我需要详细分解这些步骤，确保每一步都正确，比如在ArcGIS中使用“面转栅格”工具时，需要设置像元大小和字段，接着用“栅格转点”工具生成点要素，然后添加经纬度字段，最后导出属性表为CSV。

接下来是PostGIS的方法，引用[1]提到将各种GIS数据导入PostGIS数据库，但用户现在是要导出栅格数据为CSV。这时候可能需要使用PostGIS的函数，比如ST_PixelAsXXX系列函数来提取栅格的像素值。可能需要将栅格分解为像素点，然后导出这些点的坐标和值到CSV。但我不太确定具体的SQL语句怎么写，可能需要查阅PostGIS的文档，确认如何分解栅格并生成坐标。

另外，用户可能还需要其他工具的方法，比如GDAL的gdal2xyz.py脚本，这是一个常用的命令行工具，可以将栅格数据转换为XYZ格式，再转换为CSV。这个步骤需要用户安装GDAL，并正确使用命令行参数，比如跳过NoData值。

然后，我需要考虑用户的实际操作中可能遇到的问题。例如，在ArcGIS中添加经纬度字段时，是否正确计算了地理坐标或投影坐标？如果数据是投影坐标系，是否需要先转换为地理坐标系才能得到经纬度？或者在PostGIS中处理时，栅格的坐标系是否正确，是否需要转换？

还有，用户可能对编程方法感兴趣，比如Python使用GDAL库读取栅格数据，遍历每个像素，获取坐标和值，然后写入CSV。这需要代码示例，确保正确读取栅格的地理变换参数，计算每个像素的坐标，并处理NoData值。

在组织回答时，需要分步骤说明每个方法，确保逻辑清晰。同时，相关问题部分需要引导用户深入探讨，比如不同工具的比较、处理大型栅格的效率、坐标系的影响，以及如何处理多维栅格数据等。

需要验证每个步骤的正确性，例如在ArcGIS中导出属性表时，是否确实包含经纬度字段；使用GDAL的gdal2xyz.py时，是否正确处理了不同数据格式。可能还需要提醒用户注意数据量大的时候，CSV文件可能会很大，需要适当的内存和存储空间。

最后，确保回答符合用户的要求：结构清晰，步骤详细，引用正确，并且生成的相关问题能够覆盖用户可能进一步的需求。同时，避免使用Markdown格式，保持自然的中文表达，口语化但不啰嗦。</think>### 栅格数据转换为CSV表格的方法  
将栅格数据（如高程、温度、遥感影像等）转换为CSV文件的核心步骤是提取每个栅格像元的坐标（经纬度或投影坐标）及其对应值。以下是具体实现方法：

---

#### **1. 使用ArcGIS工具链**  
**步骤说明：**  
1. **面转栅格（可选）**  
   若需要特定区域的栅格数据，先用面状shp文件裁剪目标范围，再通过 `ArcToolbox > Conversion Tools > 转为栅格 > 面转栅格` 生成栅格数据[^2]。  
2. **栅格转点**  
   使用 `ArcToolbox > Conversion Tools > 从栅格转出 > 栅格转点`，将栅格像元转为点要素图层，每个点对应一个像元的中心坐标。  
3. **添加坐标字段**  
   在点要素属性表中，通过 `计算几何` 工具添加经度（`X`）和纬度（`Y`）字段，选择坐标系（地理坐标系则输出经纬度，投影坐标系输出平面坐标）。  
4. **导出为CSV**  
   右键点击点要素图层，选择 `导出表格`，保存为CSV格式，包含 `X`、`Y` 和栅格值字段。

---

#### **2. 使用GDAL命令行工具**  
**方法：**  
通过 `gdal2xyz.py` 脚本直接将栅格转换为文本格式（XYZ），再另存为CSV：  
```bash
gdal2xyz.py -skipnodata input_raster.tif output.xyz
# 将.xyz重命名为.csv，或使用Excel转换
```  
参数 `-skipnodata` 可跳过无效值，适用于大型栅格数据。

---

#### **3. 通过PostGIS数据库处理**  
**步骤：**  
1. **导入栅格到PostGIS**  
   使用 `raster2pgsql` 工具将栅格数据导入数据库：  
   ```bash
   raster2pgsql -s 4326 -I -C -M input.tif -F public.raster_table | psql -d your_db
   ```  
2. **提取像元坐标和值**  
   执行SQL查询，通过 `ST_PixelAsXXX` 函数分解栅格：  
   ```sql
   SELECT (ST_PixelAsPoints(rast)).*
   FROM public.raster_table;
   ```  
3. **导出查询结果为CSV**  
   使用 `COPY` 命令或PgAdmin导出工具将结果保存为CSV。

---

#### **4. Python脚本（基于GDAL库）**  
**代码示例：**  
```python
from osgeo import gdal
import csv

# 读取栅格
ds = gdal.Open("input.tif")
band = ds.GetRasterBand(1)
data = band.ReadAsArray()
gt = ds.GetGeoTransform()

# 计算坐标
with open("output.csv", "w", newline='') as csvfile:
    writer = csv.writer(csvfile)
    writer.writerow(["X", "Y", "Value"])
    for y in range(data.shape[0]):
        for x in range(data.shape[1]):
            x_coord = gt[0] + x * gt[1] + y * gt[2]
            y_coord = gt[3] + x * gt[4] + y * gt[5]
            writer.writerow([x_coord, y_coord, data[y, x]])
```  
此脚本遍历每个像元，计算其坐标并写入CSV，适用于自定义处理逻辑[^2]。

---

			
		

	



              




          




        



    





    



      

      

        
      

      





	

		
评论 4

	

  
	




  

    

      添加红包
      
    

    

      

        
        

          
          
        

        
请填写红包祝福语或标题

      

      

        
        

          
          个
        

        
红包个数最小为10个

      

      

        
        

          
          元
        

        
红包金额最低5元

      

      

        
        

          当前余额3.43元
          前往充值 >
        

      

      

        

          需支付：10.00元
        

        
        
      

    

  





  

    
    
  

  

    
    
  








  

    

      

        

          

            
            
成就一亿技术人!

          

          

        

        

          

          

            领取后你会自动成为博主和红包主的粉丝
            规则
          

        

      

      

        

          

            
              
            
            

              
hope_wisdom
 发出的红包
            

          

        

        

          

          

          

        

      

    

    

  



      
      

      
实付元

      
使用余额支付

      

        

          

            
            点击重新获取
          

        

        
扫码支付

      

      

        
          
            
          
          
        
      

      

        
        钱包余额
          0
          

            
            

              

                
抵扣说明：

                
 1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
 2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

              

            

          

      

      余额充值