Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  636 1375 (N37E16 L=160)   959    240        4 ( 556/1320- 745/1407)  -3.6 E23-E07 N35-N39
      2  158 1210 (N37E61 L=114)   739    273        0 ( 120/1175- 200/1253) -32.0 E65-E56 N34-N39
      3 1850 1020 (S18W55 L=231)   971    221        4 (1820/ 949-1885/1071)  45.0 W51-W59 S21-S15
      4  354 1001 (N21E44 L=132)   935    190        2 ( 294/ 971- 428/1056) -37.9 E48-E38 N18-N25
      5 1776  944 (S19W46 L=222)  1094    158        2 (1748/ 879-1823/1002)  24.0 W42-W51 S23-S16
      6 1431  782 (S18W17 L=193)   623    184        1 (1363/ 714-1508/ 845)  20.0 W11-W23 S20-S15
      7  197  779 (N10E66 L=110)  1463    231        4 ( 157/ 722- 224/ 833) -90.0 E71-E62 N06-N13
      8 1040  781 (S07E05 L=170)   242     56        0 (1025/ 746-1056/ 808)  56.3 E06-E03 S08-S05
      9  635  549 (S11E36 L=139)  2161    373        1 ( 577/ 383- 724/ 717)  73.1 E40-E31 S23-N00
  

    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