Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  805 1676 (N33E26 L=309)   722    213        1 ( 723/1595- 867/1733) -29.4 E34-E18 N27-N36
      2  672 1486 (N18E29 L=306)   293     91        1 ( 651/1449- 695/1530) -51.3 E31-E27 N15-N20
      3 1272 1341 (N16W11 L=346)   515    157        0 (1217/1289-1330/1403) -40.6 W06-W15 N13-N19
      4 1642 1288 (N19W37 L= 12)   335     54        0 (1611/1265-1670/1305) -11.3 W34-W39 N18-N20
      5  616  904 (S18E24 L=311)   443    119        0 ( 582/ 840- 655/ 951)  66.8 E26-E21 S21-S14
      6  214  850 (S24E60 L=275)   438     75        1 ( 184/ 818- 237/ 885) -38.7 E63-E57 S25-S21
      7  849  748 (S25E07 L=328)   513    208        1 ( 797/ 669- 886/ 847)  84.8 E10-E05 S29-S18
      8  758  522 (S41E12 L=323)   302    104        0 ( 705/ 499- 812/ 541)  -6.3 E16-E07 S41-S39
  

    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