Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  576 1474 (N43E19 L=144)   345     94        2 ( 533/1438- 613/1500)   7.1 E22-E14 N41-N44
      2  666 1231 (N25E16 L=146)   746    208        1 ( 639/1137- 689/1335)  81.3 E18-E13 N19-N32
      3  819 1238 (N22E06 L=156)   201     75        0 ( 788/1218- 860/1251)   0.0 E08-E03 N20-N22
      4  202 1050 (N25E59 L=104)   516    109        1 ( 187/ 983- 229/1097)  81.9 E60-E55 N21-N28
      5  164 1037 (N25E65 L= 98)   487    102        2 ( 148/ 995- 176/1095)  63.4 E67-E62 N22-N28
      6 1588  861 (S19W30 L=192)   269    110        2 (1545/ 821-1634/ 893)  15.9 W26-W33 S20-S18
      7 1708  649 (S37W38 L=200)   379     47        1 (1676/ 627-1730/ 673)   0.0 W34-W40 S39-S35
      8 1512  635 (S31W20 L=183)   241     94        1 (1478/ 613-1548/ 671)  18.4 W17-W23 S32-S29
      9 1250  408 (S38E04 L=159)   468    132        1 (1195/ 380-1313/ 448)   6.3 E08-W01 S39-S37
     10  545  268 (S29E64 L= 98)   743     71        2 ( 525/ 227- 559/ 307)  20.6 E74-E58 S32-S26
  

    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