Oracle和Postgis資料庫地理座標系下面積計算
阿新 • • 發佈：2021-07-21
要求：計算投影面積，單位公頃，保留4位小數
Oracle
round(
(st_area(
st_transform(shape,
(
floor((st_x(st_centroid(shape))+1.5)/3)
-25+4513))
)
/10000),
4)
Postgis
round(
(st_area(
st_transform(geom,
(
floor((st_x(st_centroid(geom))+1.5)/3)
-25+4513)::int4)
)
/10000)::numeric,
4)
計算思路：
1、獲取圖形中心點X座標
Centroid_X = st_x(st_centroid(geom))
2、計算3度帶帶號
3_degree=floor((Centroid_X+1.5)/3)
3、計算3度帶帶號對應SRID
| SRID | COORD_REF_SYS_NAME | X |
|---|---|---|
| 4513 | CGCS2000 / 3-degree Gauss-Kruger zone 25 | 75 |
| 4514 | CGCS2000 / 3-degree Gauss-Kruger zone 26 | 78 |
| 4515 | CGCS2000 / 3-degree Gauss-Kruger zone 27 | 81 |
| 4516 | CGCS2000 / 3-degree Gauss-Kruger zone 28 | 84 |
| 4517 | CGCS2000 / 3-degree Gauss-Kruger zone 29 | 87 |
| 4518 | CGCS2000 / 3-degree Gauss-Kruger zone 30 | 90 |
| 4519 | CGCS2000 / 3-degree Gauss-Kruger zone 31 | 93 |
| 4520 | CGCS2000 / 3-degree Gauss-Kruger zone 32 | 96 |
| 4521 | CGCS2000 / 3-degree Gauss-Kruger zone 33 | 99 |
| 4522 | CGCS2000 / 3-degree Gauss-Kruger zone 34 | 102 |
| 4523 | CGCS2000 / 3-degree Gauss-Kruger zone 35 | 105 |
| 4524 | CGCS2000 / 3-degree Gauss-Kruger zone 36 | 108 |
| 4525 | CGCS2000 / 3-degree Gauss-Kruger zone 37 | 111 |
| 4526 | CGCS2000 / 3-degree Gauss-Kruger zone 38 | 114 |
| 4527 | CGCS2000 / 3-degree Gauss-Kruger zone 39 | 117 |
| 4528 | CGCS2000 / 3-degree Gauss-Kruger zone 40 | 120 |
| 4529 | CGCS2000 / 3-degree Gauss-Kruger zone 41 | 123 |
| 4530 | CGCS2000 / 3-degree Gauss-Kruger zone 42 | 126 |
| 4531 | CGCS2000 / 3-degree Gauss-Kruger zone 43 | 129 |
| 4532 | CGCS2000 / 3-degree Gauss-Kruger zone 44 | 132 |
| 4533 | CGCS2000 / 3-degree Gauss-Kruger zone 45 | 135 |
srid=3_degree-25+4513
4、地理座標系轉投影座標系
Oracle
geometry=st_transform(shape,srid)
Postgis
geometry=st_transform(shape,srid::int4)
注：srid需要型別轉換，否則報錯誤: 函式 st_transform(geometry, double precision) 不存在
因為floor()結果為double型別，postgis支援的函式為st_transform(geometry, int4)
5、計算面積
area=st_area(geometry)
6、平方米轉公頃，保留4位小數
Oracle
area_hectares=round(area/10000,4)
Postgis
area_hectares=round((area/10000)::numeric,4)
注：(area/10000)需要型別轉換，否則報錯誤: 函式 round(double precision, integer) 不存在
因為(area/10000)為double型別，postgis支援的函式為round(numeric, int4)