2017-02-06 10 views
2

私は約100,000の緯度と経度のペア(4小数点以下まで)を持っています。私はそれぞれのペアを米国の州に割り当てたいと思います。誰でもSASでこれを行う方法を知っていますか?このタスクのためにシェープファイルをSASにインポートすることは可能ですか?SASで緯度と経度の状態をマップする

+0

SASにはシェープファイルが事前にロードされていますが、SASの特定の投影にlat/longを投影する必要があります。地図投影法を使い、それを操作するのはどれくらい快適ですか? – Joe

+0

SASで地図投影法を使ったことはありません。私はSAS Baseの手順でかなり経験がありますが、プリロードされたシェイプファイルではありません。 –

+0

Hmm。私はそれを見て、 "投影SAS gmap gproject" - それがあなたに合っているかどうかを見てみましょう。私はあなたに実際の答えを与えるのに十分よく知っているとは思わないが、それはあなたが行く必要のある方向だと思う。 – Joe

答えて

3

注:以下の回答は、SAS Graphのライセンスがあり、インストール時にマップライブラリが適切に設定されていることを前提としています。 SAS University Editionでは動作しません。 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; 
+1

Reezaありがとうございます。これは私が必要としていたものです。私はジンサイドについて知らなかった。ありがとうございました! –

関連する問題