Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1509 1733 (N30W58 L=202)  2508    448        5 (1402/1668-1610/1768) -10.8 W47-W68 N24-N32
      2 1033 1051 (N06W01 L=145)   333    104        0 ( 995/ 996-1085/1098)  23.2 E02-W05 N03-N07
      3 1183 1029 (N01W09 L=153)   274     90        0 (1139/ 988-1225/1055) -53.1 W06-W10 S02-N02
      4  320  986 (N18E46 L= 98)   327     70        1 ( 306/ 951- 336/1034) -90.0 E47-E44 N16-N21
      5 1836  895 (S26W51 L=195)   717    141        1 (1800/ 822-1867/ 960)  41.2 W46-W55 S29-S22
      6 1580  701 (S29W26 L=171)   598    156        0 (1526/ 619-1638/ 784)  38.7 W21-W31 S33-S25
      7 1859  749 (S36W53 L=197)   482     97        0 (1826/ 726-1883/ 774)   9.5 W49-W56 S37-S34
      8 1288  696 (S21W07 L=151)   289     85        0 (1246/ 675-1341/ 711)  -8.1 W03-W10 S21-S19
      9 1161  500 (S30E06 L=138)   345     85        0 (1119/ 469-1192/ 536)  18.4 E09-E03 S30-S27
     10  564  435 (S18E46 L= 98)   288     19        0 ( 544/ 414- 584/ 453) -63.4 E47-E44 S20-S17
     11  645  342 (S26E46 L= 98)   717    190        1 ( 618/ 276- 663/ 405)  52.1 E50-E42 S31-S22
     12  501  306 (S26E60 L= 84)  1578    272        1 ( 475/ 239- 517/ 372)  42.3 E66-E54 S30-S20
  

    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