Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1281 1805 (N52W15 L=261)   750     70        3 (1255/1760-1305/1841)  76.0 W11-W17 N47-N57
      2  558 1659 (N32E43 L=204)   936    122        4 ( 528/1581- 581/1714) -55.0 E46-E39 N26-N36
      3 1484 1480 (N28W27 L=274)   445    119        2 (1440/1442-1534/1503)  23.2 W24-W31 N25-N29
      4  269 1380 (N12E58 L=188)   837    125        2 ( 229/1321- 299/1432) -30.3 E65-E53 N08-N15
      5 1734 1179 (N12W46 L=293)   527     59        0 (1698/1146-1759/1206) -21.8 W43-W48 N09-N13
      6 1691 1052 (N03W43 L=290)   281     37        0 (1678/1019-1705/1077) -71.6 W42-W44 N01-N04
      7 1710  640 (S20W54 L=301)   512     66        2 (1685/ 607-1737/ 669)  63.4 W52-W56 S22-S18
      8 1357  645 (S25W26 L=272)   243     46        1 (1335/ 622-1382/ 665)  45.0 W24-W27 S26-S23
      9 1234  552 (S32W20 L=266)   224     45        1 (1216/ 531-1256/ 571)  71.6 W18-W21 S33-S30
     10 1542  358 (S41W56 L=303)  1077     87        1 (1503/ 315-1568/ 413) -35.0 W50-W64 S43-S36
     11  575  309 (S57E42 L=205)  1269    128        3 ( 522/ 275- 632/ 349)  16.9 E55-E32 S60-S53
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y