Lunarpedia talk:Autostub1 test 1

Erratum

• Category tag for 'Location' should have been 'Locations'

-- Strangelv 18:27, 4 November 2006 (PST)

Problems

• How to tie the satellite crater articles to the main crater articles when the USGS database doesn't link them?

• What is the area of each sector in km**2?

• How might the script determine if a given feature is in a given sector? If we treat each item as a perfect circle, how do we determine which sectors include portions of that circle?

• Can we automatically do maps for features near the poles?

-- Strangelv 18:13, 4 November 2006 (PST)

To Do

• Sector areas ---- I posted the solution to this Charles F. Radley 17:53, 16 January 2007 (PST)
• Match sectors and features (not just centers) -- easy to solve, will do on a weekend or snow day Charles F. Radley 17:53, 16 January 2007 (PST)
• Match satellite craters with their central craters - Where is the list of satellite craters? Charles F. Radley 17:53, 16 January 2007 (PST)

I have the data downloaded, but I can't find the URL from the USGS page I got it from. -- Strangelv 18:38, 16 January 2007 (PST)

• Border case handling for when feature center is on the boundary of two or four sectors - What do you want it to do in those cases ? Charles F. Radley 17:53, 16 January 2007 (PST)

Possibly list as included in all of the above? -- Strangelv 18:38, 16 January 2007 (PST)

• Automatically load and parse tab delimited CSV --- Is this done now ? Charles F. Radley 17:53, 16 January 2007 (PST)

No. You may notice the hard-coded data entry in the script. 8) I basically cranked out the essentials to produce the test page and then got sidetracked as well as stopped by a few technical issues listed on this page. -- Strangelv 18:38, 16 January 2007 (PST)

• Automatically output articles in XML for bulk import
• Automatically generate map files for suitably large features
• Automatically generate sector articles and maps

-- Strangelv 18:13, 4 November 2006 (PST)

Bonus Points

• Is there a public domain database of rille?
• Are there other public domain databases of use that we might be able to merge into this to add details to sectors and features (such as surface samples)?

-- Strangelv 18:13, 4 November 2006 (PST)

The precise coordinates of the Apollo sites are available on the web link provided on that page. Charles F. Radley 17:50, 16 January 2007 (PST)

Just need to add those to the database if they're not already there. I believe several features named on the ground are already in there. -- Strangelv 18:38, 16 January 2007 (PST)

• Is there a public domain database that represents non-circular features in a non-circular manner?

-- Strangelv 16:25, 27 December 2006 (PST)

How about approximating them to ellipses? Standard geometry equations. Charles F. Radley 17:49, 16 January 2007 (PST)

Sounds great if it can be found. I've only seen one database though, and it treats everything as circles whether perfectly appropriate or wildly inappropriate. Manually changing them ourselves at this stage isn't feasible given the sheer number of features in the database. -- Strangelv 18:38, 16 January 2007 (PST)

Software Development Plan

Thanks for the responses to my various questions. Essentially what I would like to see here is a detailed requirements specification and/or a detailed design document. that will address your valid concerns that "we need to get it right first time"...

In my experience, the most common and most devastating problem encountered by software development projects is.... screwed up requirements, ie requirements which are vague, ambiguous, confusing, misleading, impossible, contradictory, meaningless, non-sensical, and/or incomplete, or just plain wrong. So that is why I am asking lots of questions to make sure the requirements are crystal clear, and we know exactly what we are trying to achieve before we go zooming up the wrong road.

Thanks.

Charles F. Radley 19:41, 16 January 2007 (PST)

Thank you for pointing out this blazingly obvious point that mysteriously never occurred to me. In my defense, I've been working on an iterative development concept (at least until I got detoured) that is compatible with cranking out something now and getting a clue later (even if it means scrapping or shelving a lot of code in the process).

How confusing is my internal documentation? A historical complaint is that my commenting has often been less than fully helpful.

This probably requires its own page, as we're moving beyond the first test. How about Lunarpedia:Autostub1 (...which already exists -- I forgot all about it)?

-- Strangelv 01:41, 17 January 2007 (PST)

Sector areas code:Set up variables to generate article stub Are we ready to incorporate the equations for sector area into the script? Charles F. Radley

Sector Areas could be precomputed and referenced by latitude, as there's basically only five or so different areas unless we want to get really serious and go by actual complex geometry instead of treating it as a sphere (or for partial bonus credit, a slightly flattened sphere) (for example, I've been told that Earth is slightly pear shaped).

Presumably the computation or selection of area could be added to 'Set up variables to generate article stub' and then added to the template graphics below that before to_out is run through tidypipes(), as required by {{Circular Feature}}.

--Strangelv 01:41, 17 January 2007 (PST)

James, what do you mean by "five or so different areas " ? Do you mean there are only 5 sectors?

The solution I posted assumes the Moon is a perfect sphere, for our purposes there is no reason to use more ellaborate models. I believe the equations I provided are adequate. We are not trying to perform a pinpoint landing of a spacecraft, nor make topographic measurements, we are only giving the users some visual experience.

Please provide the lattitude/longitude boundaries of each sector and I will calculate the area by hand. Where do you want me to put the results?

Charles F. Radley 05:02, 17 January 2007 (PST)

24 sectors along the equator have the same area; 2 polar sectors have the same area; all other sectors share the same area with 47 other sectors. Sorry if 'area' is a confusing term to use here. -- Strangelv 06:21, 17 January 2007 (PST)

Source Code

Warning: this was getting mauled in a few places. I've tinkered with the formatting to make it work.

And yes, the code is really rough in a few places (like where reinventing the wheel is faster than remembering where the wheel is). The language is Python. 2.4 and later should work; earlier may require replacing my '//' division with integer coercion.

 #%#%#%#%#%#%#%#%#%#%#%#%#%#%#%#
 #                             #
 # Feature Autostub Generator  #
 #                             #
 # Licence pending             #
 #                             #
 #%#%#%#%#%#%#%#%#%#%#%#%#%#%#%#

 def dostub(thisentry):
     """Create stub article from database
     """
     #   
     # Set up variables to generate article stub
     #
     #%#%#%#%#%#%#%#%#%#%#%#%#%#%%#%#%#%#%#
                                          #
     name = thisentry[2]                  # name of feature
     mapname = makemapname(name, None)    # filename for this feature's
                                          #      automatically generated
                                          #      map from Clementine data
     Lat = float(thisentry[3])            # USGS latitude
     lat = do_lat(Lat)                    # string with proper latitude
     me_long = poslong(thisentry[4])      # longitude (Mean Earth)
     mo_long = me2mo(me_long)             # longitude (Mare Orientale)
     me_coord = coordstr(me_long, lat)    # Mean Earth coordinate string
     mo_coord = coordstr(mo_long, lat)    # Mare Orientale coordinate string
     origin = thisentry[24]               # 'Origin' column in database
                                          #          warning: this seems to
                                          #          be almost anything
     featuretype = thisentry[23]          # Is this a mare, a crater, other?
     dia = thisentry[10]                  # feature diameter, if applicable
     approvalStat = thisentry[18]         # year approval information
     approvalYear = thisentry[19]         # year of approval (if approved)
     reference = thisentry[21]            # 'reference' column in database
     inSector = isinsector(mo_long, Lat)  # what sector is this feature
                                          #      centered in?
     sectors = sectorlist(mo_long, lat, dia)
                                          # -----what sectors include it?
     typecat = dotypecat(featuretype)     # simplified category type
     domap(mapname, mo_long, lat, dia)    # create the map
     localeDescr = descrlocale(mo_long, Lat, typecat)
                                          # -----text description of location
     approvalStr = doapproval(approvalStat, approvalYear)
                                          # -----text description of approval,
                                          #      if approved.
     #%#%#%#%#%#%#%#%#%#%#%#%#%#%%#%#%#%#%#
     #
     # start generating the article here:
     #
     to_out  = []
     to_out += ["{{Circular Feature "]
     to_out += ["name        = "+name]
     to_out += ["featuretype = "+typecat]
     to_out += ["origin      = "+origin]
     to_out += ["map         = [[Image:"+mapname+"|240px]]"]
     to_out += ["mo_coord    = "+mo_coord]
     to_out += ["me_coord    = "+me_coord]
     to_out += ["dia         = "+dia]
     to_out += ["approval    = "+approvalStr]
     
     # Tidily install pipes for the table 
     to_out = tidypipes(to_out)
     
     to_out += ["}}"]
     to_out += ["'''"+name+"''' "+localeDescr]
     to_out += ["<BR/><BR/>\n\n"]
     # if inSector !="AAAAAAAAA!":
     to_out += ["The center of "+name+" is in sector [[Sector_"+inSector+"]]"]
     # else:
     #      inSectors = something()
     to_out += ["<BR/><BR/>\n\n"]
     to_out += ["==See Also=="]

     for q in sectors:
         to_out += ["[["+q+"]]"]
     
     to_out += ["<BR/><BR/>\n\n"]
     to_out += ["==References=="]
     to_out += [refbreak(reference)]
     to_out += ["{{Autostub}}"]
     to_out += ["[[Category:Location]]"]
     to_out += ["[[Category:"+typecat+"]]"]
     to_out += [""]

     return to_out

 # def main():
 #     # load database
 #     didxml = templatetop()
 #     for q in database: didxml += dostub(q) + templatemid(); print ".",
 #     didxml += templatefin()
 #     # save didxml

 def makemapname(name, other):
     """
     """
     name2=""
     for q in name:
         if q == " ": name2 += "_"
         else: name2 += q
     if other: name2 += "_"+other
     name2 += ".png"
     return name2

 def poslong(long_):
     """convert string of possibly noncompliant Mean Earth longitude into float or fixed of compliant Mean Earth longitude
     """
     long_ = float(long_)
     if long_ < 0 : long_ += 360
     return long_

 def me2mo(me_long):
     """
     """
     mo_long=float(me_long)-267.2
     if mo_long <0: mo_long += 360.0
     #print "me_long == "+str(me_long)
     return mo_long

 def coordstr(long_,lat):
     """ temporary implementation
     """
     return str(lat)+" "+str(long_)

 def descrlocale(long_, lat, what):
     """quickie hack
     """
     if   lat > 5: vert = "a northern"
     elif lat < 5: vert = "a southern"
     else: vert = "an equatorial"

     if   long_ <= 165 and long_ >= 15: side = "near side"
     elif long_ <= 195 and long_ >= 345: side = "far side"
     elif long_ >= 165 and long_ <= 195: side = "east edge"
     else: side = "west edge"
     
     return "is "+vert+" "+what+" on the "+side+" of Luna."

 def isinsector(long_,lat):
     """can we determinewhat sector you're in?
     """
     # special case handling?
     def AAAAAAAAA():
         # border case handling
         return "AAAAAAAAA!"
     # find latitude part or bust
     if lat == 0.0: LT="E" # this is equatorial
     else:
          if lat > 0: LTH = "N"
          else:
               lat = -lat
               LTH = "S"
          if ((lat+7.5)//15) == ((lat+7.5)/15.0): return AAAAAAAAA()
          LT = str(int(((lat+7.5)//15)*15)) + LTH
          if LT == "90": return LTH # this is a polar sector
     # find longitude part or bust
     if long_ < 7.5 or long_ > 352.5: LG = "000" # prime meridian sector
     else:
          if ((long_+7.5)//15) == ((long_+7.5)/15.0): return AAAAAAAAA()
          LG = str(int(((long_+7.5)//15)*15))
          if len(LG)<3: LG = "0"+LG
     # return the result!
     return LT+LG

 def dotypecat(stringie):
     """
     """
     if stringie == "Rima, rimae": return "rima"
     # Most definitely need to list the complete set
     return stringie

 def sectorlist(long_,lat,dia):
     """
     """
     return ["sectorlist() not implemented"]

 def domap(mapname,long_,lat,dia):
     """
     """
     return "domap() not implemented"

 def doapproval(approvalStat, approvalYear):
     """temporary hack; needs to acconut for unapproved
     """
     return approvalStat+" in "+approvalYear

 def tidypipes(table):
     """takes a list oy strings, finds the longest one, and tidily adds pipes to the right of each line.

     maxlen: maximum length
     table:  input list of strings
     table2: output list of strings
     """
     maxlen = 0
     for q in table:
         if len(q)>maxlen: maxlen = len(q)
     table2 = []
     for q in table:
         qq = q
         while (len(qq)<maxlen):
               qq += " "
         table2 += [qq+" |"]
     return table2

 def do_lat(lat):
     """returns string for latitude
     """
     if lat > 0: return str(lat)+"N"
     if lat == 0: return "E"
     return str(-lat)+"S"

 def replacer(stringie, old, new):
     """
     """
     name2=""
     for q in stringie:
         if q == old: name2 += new
         else: name2 += q
     return name2

 def refbreak(ref):
     """
     """
     return replacer(ref,";",";<BR/>")

 ##def ():
 ##    """
 ##    """
 ##    return "() not implemented"
 ##
 ##def ():
 ##    """
 ##    """
 ##    return "() not implemented"

 #Item = ["Earth","Moon","Rima Diophantus","29.0","-33.0"," "," "," "," ","east","150.0","EU","Europe","GR","Greek","LTO39B2","LTO-39B2","5","Adopted by IAU","1985","0","null","RI","Rima, rimae","Named from nearby crater."]
 Item = ["Earth","Moon","Rima Cardanus","11.4","-71.5"," "," "," "," ","east","175.0","EU","Europe","IT","Italy","LOC-1","LOC-1","5","Adopted by IAU","1964","67","The System of Lunar Craters, Quadrants I, II, III, IV; by D. W. G. Arthur and others; Communications of the Lunar and Planetary Laboratory, vol. 2, no. 30, 1963; vol. 3, no. 40, 1964; vol. 3, no. 50, 1965; vol. 5, no. 70, 1966.","RI","Rima, rimae","Named from nearby crater."]

 What=dostub(Item)
 for q in What: print q