Широта и долгота в SAS

Question

У меня есть около 100 000 пар широты и долготы (до 4 десятичных знаков), и я хотел бы назначить каждую пару в состояние США. Кто-нибудь знает, как это сделать в SAS? Можно ли импортировать шейп-файлы в SAS для этой задачи?

Answer 1

Примечание: следующий ответ предполагает, что у вас есть лицензия на SAS Graph, и библиотека Карт правильно настроена при установке. Это не будет работать в SAS University Edition. Если вам необходимо загрузить файлы MAPS, они находятся здесь: http://support.sas.com/rnd/datavisualization/mapsonline/index.html

Вам не нужно импортировать файл SHP для США, SAS уже имеет эти встроенные устройства. Вы можете использовать PROC GINSIDE, чтобы определить, какое состояние и/или округ, точки расположены.

Пример расположен здесь: https://support.sas.com/documentation/cdl/en/grmapref/69722/HTML/default/viewer.htm#p0qjcc8hugcjb2n1x3bmuaar16f0.htm

И скопировал здесь, для SO правил.

goptions reset=global border; 
data gpscounties; 
    input longitude latitude site $; 
    x=longitude*arcos(-1)/180; 
    x=x*(-1); 
    y=latitude*arcos(-1)/180; 
datalines; 
-77.0348 40.0454 a 
-78.4437 39.1623 b 
-78.4115 39.3751 c 
-78.7646 40.6354 d 
; 
run; 
proc ginside data=gpscounties map=mapssas.counties out=gpscounties; 
    id state county; 
run; 
proc sort data=gpscounties; 
    by site; 
run; 

proc print data=gpscounties; 
    var site state county x y; 
run; 
quit;