Я получил файл формы из http://www.savgis.org/Ecuador/provincias.zip без CRS и хочу добавить его. Я обнаружил, что для моих данных это значения EPSG: 24878, PSAD56 / UTM zone 18S https://epsg.io/24878. Я сделал это
Ecuador<- read_sf("provincias.shp")
st_crs(Ecuador) # no CRS
Ecuador=st_set_crs(Ecuador,24878) # add crs to the data
st_crs(Ecuador) # works
st_area(Ecuador) # area is not right
Я проверил, что геометрия <MULTIPOLYGON [m]>
, я хочу знать, как ее изменить или как мне добавить правильные блоки CRS. Конечная область должна быть около 250747.1 [km^2]
Спасибо