Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  309 1527 (N34E67 L=193)   796     81        1 ( 289/1495- 328/1581) -26.6 E73-E62 N32-N37
      2 1116 1444 (N30W07 L=267)   375    173        0 (1038/1417-1190/1481) -21.8 W01-W11 N28-N32
      3  199 1184 (N12E64 L=196)  1139    122        4 ( 160/1138- 228/1243)  54.5 E69-E60 N08-N15
      4 1889 1037 (N02W68 L=328)  3519    293       14 (1820/ 917-1943/1119) -38.7 W58-W79 S05-N07
      5  750 1001 (N02E17 L=243)   204     81        1 ( 712/ 972- 775/1026) -36.9 E19-E15 N00-N03
      6  235  923 (S04E58 L=203)   646    115        0 ( 216/ 867- 259/ 973)  54.5 E60-E54 S07-N00
      7  367  805 (S11E46 L=215)   637    192        0 ( 313/ 739- 396/ 870)  56.3 E51-E42 S15-S06
      8  183  733 (S17E70 L=190)   883    106        2 ( 156/ 688- 213/ 774)  22.4 E81-E63 S20-S13
      9  445  517 (S30E46 L=214)  2696    451        4 ( 397/ 361- 515/ 671)  39.6 E65-E35 S43-S19
     10 1724  549 (S29W58 L=319)   772    148        0 (1695/ 488-1748/ 613) -48.4 W56-W64 S33-S24
     11 1475  519 (S30W33 L=294)  2221    466        3 (1238/ 451-1614/ 605)  13.6 W15-W45 S34-S24
  

    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