Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1265 1868 (N59E20 L=271)   624     45        1 (1246/1836-1281/1891) -31.0 E25-E15 N55-N63
      2 1098 1483 (N21E09 L=281)   273    120        0 (1061/1426-1122/1519) -29.7 E12-E05 N18-N22
      3 1789 1129 (N23W47 L=337)  1745    357        2 (1721/1035-1860/1241) -28.1 W38-W56 N17-N27
      4 1578 1147 (N16W29 L=319)   337    113        1 (1541/1106-1602/1199)  76.0 W27-W30 N12-N20
      5 1651 1001 (N11W38 L=328)   525     86        1 (1633/ 945-1669/1049) -63.4 W35-W39 N07-N13
      6  247  962 (S29E54 L=237)  1797    272        2 ( 131/ 938- 390/ 995)   5.2 E72-E39 S31-S26
      7 1832  748 (N04W64 L=354)   748    104        0 (1775/ 729-1882/ 765)  19.7 W57-W71 N01-N06
      8  862  710 (S29E00 L=290)   382     97        1 ( 818/ 677- 898/ 749)  71.6 E02-W01 S31-S25
      9 1715  648 (S05W57 L=347)   323     28        1 (1696/ 625-1729/ 672) -71.6 W54-W58 S06-S03
     10 1302  544 (S24W33 L=323)   646    227        0 (1199/ 521-1408/ 565)  33.7 W26-W38 S27-S19
     11 1557  418 (S22W59 L=349)  1871    230        6 (1497/ 340-1625/ 502) -63.4 W54-W64 S27-S15
  

    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