Zakładając, że twoje obszary zainteresowania są raczej małe w porównaniu do globu, możesz ustawić niestandardową projekcję poprzecznego mercatora.
Musisz znać współrzędne geograficzne lat_0 i lon_0 początku twojego CRS, a także kierunek osi xiy:
+proj=tmerc +lat_0=51.4 +lon_0=7 +k=1 +x_0=0 +y_0=0 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs
Do użytku z Arcgis plik .prj powinien zawierać:
PROJCS["Transverse_Mercator",GEOGCS["GCS_WGS_1984",DATUM["D_unknown",SPHEROID["WGS84",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",51.4],PARAMETER["central_meridian",7],PARAMETER["scale_factor",1],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1]]
Jeśli nie są one równoległe do kierunków północnego i wschodniego, konieczny byłby obrót za pomocą skośnego rzutu mercatora:
+proj=omerc +lat_0=51.4 +lonc=7 +alpha=-10 +k=1 +x_0=0 +y_0=0 +gamma=0 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs
Plik .prj dla tej projekcji to:
PROJCS["Hotine_Oblique_Mercator_Azimuth_Center",GEOGCS["GCS_WGS_1984",DATUM["D_unknown",SPHEROID["WGS84",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Hotine_Oblique_Mercator_Azimuth_Center"],PARAMETER["latitude_of_center",51.4],PARAMETER["longitude_of_center",7],PARAMETER["azimuth",-10],PARAMETER["scale_factor",1],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1]]
Oto oba razem (tmerc na niebiesko i omerc na czerwono):
Lat_0 i lon_0 (lub lonc) są współrzędnymi lokalnego pochodzenia CRS. Musisz je eksplorować za pomocą GPS lub z innych map. Kąt, który musisz obliczyć na podstawie współrzędnych dwóch znanych punktów na lokalnej osi X.
Jedyne, co nie działa zgodnie z oczekiwaniami, to milimetr jako jednostki. Powinien działać, ale nie działa w obecnym QGIS :-(
http://hub.qgis.org/issues/9414
Aktualizacja
Obecne wersje QGIS są teraz w stanie pracować z lokalnymi rzutami przy użyciu mm jako jednostek.