code: def sample_lonlat(self, n):
          
          # From http://en.wikipedia.org/wiki/Ellipse#General_parametric_form
          # However, Martin et al. (2009) use PA theta "from North to East"
          # Definition of phi (position angle) is offset by pi/4
          # Definition of t (eccentric anamoly) remains the same (x,y-frame usual)
          # In the end, everything is trouble because we use glon, glat...
  
          radius = self.sample_radius(n)
          a = radius; b = self.jacobian * radius
  
          t = 2. * np.pi * np.random.rand(n)
          cost,sint = np.cos(t),np.sin(t)
          phi = np.pi/2. - np.deg2rad(self.theta)
          cosphi,sinphi = np.cos(phi),np.sin(phi)
          x = a*cost*cosphi - b*sint*sinphi
          y = a*cost*sinphi + b*sint*cosphi
          
          if self.projector  is None:
              logger.debug("Creating AITOFF projector for sampling")
              projector = Projector(self.lon,self.lat,'ait')
          else:
              projector = self.projector
          lon, lat = projector.imageToSphere(x, y)
          return lon, lat
docs: Sample 2D distribution of points in lon, lat