I am trying to convert latitude and longitude points to x/y coordinates, generate random variates from a bivariate normal distribution, and then convert those random x/y points back into latitude longitude coordinates. I am trying to get them to fall between 3 standard deviations of the center point that I pick, but I don't think that it is working. Here is my code so far. 

```
    XYPos = GeoGridPosition[GeoPosition[{46.6572, -41.3552}], "Mercator"]
    mu1mu2 = Extract[XYPos, {1}]
    RandomXYs = 
      RandomVariate[BinormalDistribution[mu1mu2, {30/3, 30/3}, 0], 1000]
    RandomLatLongs = GeoPosition[GeoGridPosition[RandomXYs, "Mercator"]]

```

I am trying to get the 1000 points that I generate to fall within a circle 30 nautical miles in radius, clearly my units are off. Any help would be appreciated!