Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1052 1858 (N57E29 L= 35)  1921    196        3 ( 947/1826-1145/1882)  29.4 E38-E20 N51-N62
      2  838 1508 (N23E23 L= 42)   230     93        0 ( 811/1486- 876/1528)   0.0 E25-E19 N22-N24
      3 1497 1207 (N19W26 L= 90)   977    282        5 (1359/1158-1591/1252)   3.2 W15-W33 N15-N22
      4 1647 1199 (N22W38 L=102)   339     64        1 (1622/1169-1678/1240)  78.7 W35-W40 N19-N24
      5 1643 1106 (N16W38 L=103)   769    161        4 (1580/1037-1711/1153)  48.0 W34-W44 N10-N20
      6 1042  693 (S21W08 L= 73)   230    105        2 (1000/ 666-1078/ 723)  59.0 W05-W09 S23-S18
      7 1709  649 (S09W56 L=121)   466     54        1 (1679/ 627-1738/ 679)  36.9 W54-W59 S10-S06
      8 1708  575 (S13W61 L=125)   578     71        0 (1683/ 540-1729/ 601) -38.7 W57-W64 S15-S11
      9 1364  500 (S26W36 L=100)   905    209        1 (1288/ 469-1460/ 555)  41.6 W30-W40 S28-S20
     10 1063  270 (S50W28 L= 93)   310     42        0 (1048/ 249-1085/ 293) -21.8 W25-W30 S51-S48
     11 1264  219 (S48W52 L=116)  1290    150        2 (1205/ 188-1309/ 263)   4.1 W44-W61 S50-S43
  

    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