Sorry if my way of writing is wrong, I'm a data scientist and I have little to none experience in GIS.

I would like to use an open dataset (from http://opendatastore.brussels/fr/ to be precise). It contains mobility data and some locations are encoded in EPSG:31370 (Belgian Lambert 72). I'd like to transform these locations in GPS coordinates (EPSG:4326 if I'm not mistaken).

I'm not sure to understand how you could derive a formula to automatically transform data from EPSG:31370 to EPSG:4326.

Does it exist Python libraries (or R or any other languages) in order to automate these location transformations? I found an interesting website (https://epsg.io/transform) to do it manually but I have +/- 250.000 transformations to do.

  • 2
    What is the intended output? A shapefile? A gpx-track? An excel-list? – Erik Dec 5 '19 at 16:01
  • An excel list. I already have points in GPS coordinates and I'd like to incorporate all of them in a single list. Maybe at the end, I'd like to print all these points on a map of Belgium. – Dust009 Dec 5 '19 at 16:18
  • You could download QGIS, open the dataset in QGIS, then reproject layer to 4326, then export to csv or xlsx. – jbalk Dec 5 '19 at 18:09
  • I have a Python solution in this answer here: gis.stackexchange.com/a/339783/122597 You will have to change this line: Transformer.from_crs(crs_from=31370, crs_to=4326) – Stefan Dec 5 '19 at 18:34

If I understand your excellent question, the solution is quite simple. 1) Starting with a blank QGIS project, add the Belgian Lambert 72 layer to QGIS. 2) In the Table of Contents, output a new 4326 layer by right-clicking the Belgian Lambert 72 layer and selecting Save As... 3) the resulting window will provide the opportunity to choose the output layer's CRS (in your case 4326). Voila!, the output layer is 4326!

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.