Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  527 1718 (N47E51 L=138)  1365    325        6 ( 484/1666- 561/1762) -26.6 E59-E45 N42-N50
      2  314 1390 (N22E56 L=134)  1053    227        3 ( 298/1332- 343/1480) -56.3 E60-E51 N17-N27
      3  394 1256 (N14E45 L=145)   695    219        4 ( 353/1194- 456/1315)  66.8 E47-E40 N09-N17
      4  261 1275 (N14E58 L=132)  1857    386        9 ( 183/1221- 329/1311)   6.7 E68-E51 N10-N16
      5  231 1196 (N09E60 L=130)  1114    256        3 ( 195/1162- 295/1252)  -8.7 E65-E52 N06-N12
      6  247  854 (S12E56 L=133)  1706    242        4 ( 204/ 791- 299/ 926) -41.2 E61-E51 S15-S07
      7 1762  832 (S07W54 L=244)   410     75        1 (1741/ 799-1778/ 864) -53.1 W51-W55 S09-S05
      8  888  823 (S10E08 L=182)   216     83        0 ( 849/ 791- 924/ 860)  45.0 E10-E05 S12-S07
      9 1099  450 (S35W08 L=198)   530     91        2 (1051/ 430-1150/ 474)   7.1 W04-W12 S36-S33
     10  722  418 (S39E21 L=168)  1086    289        4 ( 601/ 379- 835/ 454)  15.3 E33-E11 S43-S35
  

    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