Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  370 1533 (N30E52 L=269)  3122    347        5 ( 312/1382- 422/1663) -30.1 E83-E44 N19-N41
      2 1440 1290 (N13W26 L=347)   264     90        0 (1411/1260-1475/1327) -45.0 W24-W28 N11-N15
      3 1411 1135 (N04W24 L=344)   563    143        2 (1377/1082-1452/1201)  66.8 W21-W26 N00-N07
      4  748  904 (S10E17 L=304)   209     72        0 ( 715/ 887- 788/ 920) -21.8 E19-E14 S11-S09
      5  398  835 (S14E42 L=278)   460    158        1 ( 372/ 766- 450/ 877)  45.0 E45-E37 S18-S11
      6  766  821 (S16E16 L=305)   320     65        1 ( 730/ 796- 797/ 844) -26.6 E18-E14 S17-S14
      7 1168  771 (S19W09 L=330)   188     80        1 (1152/ 726-1194/ 803)  90.0 W08-W10 S21-S16
      8 1714  649 (S25W53 L= 13)   310     50        1 (1702/ 624-1732/ 674) -31.0 W50-W55 S26-S23
      9 1108  446 (S40W07 L=327)   260     80        0 (1068/ 431-1149/ 457)   9.5 W03-W09 S41-S39
     10  936  227 (S60E10 L=310)  1544    217        4 ( 840/ 190-1038/ 267)  11.8 E23-W01 S63-S55
  

    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