Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1755 1463 (N09W63 L=311)  2307    363        8 (1675/1415-1813/1532) -38.7 W55-W71 N05-N13
      2 1756 1257 (S02W53 L=302)   379     88        1 (1729/1217-1790/1296) -90.0 W51-W55 S04-N00
      3  121 1168 (N31E64 L=185)   875    164        1 (  99/1135- 149/1194) -63.4 E66-E60 N29-N33
      4  685 1009 (N10E20 L=229)   286     88        2 ( 639/ 985- 725/1023)   0.0 E22-E17 N08-N11
      5  861  896 (N00E12 L=237)   825    170        2 ( 770/ 869- 938/ 935) -24.0 E16-E07 S02-N04
      6  195  614 (S03E75 L=173)  1151    196        2 ( 164/ 592- 228/ 636) -15.3 E83-E70 S04-S01
      7  482  511 (S14E48 L=200)   296     55        2 ( 462/ 490- 508/ 536) -56.3 E49-E46 S15-S12
      8  546  458 (S18E46 L=203)   663    180        4 ( 511/ 413- 587/ 514)  32.0 E49-E41 S21-S15
      9  318  456 (S14E68 L=181)  2240    255        3 ( 286/ 388- 348/ 500)  81.9 E74-E63 S19-S11
  

    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