Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  679 1729 (N26E46 L= 89)  1590    624        5 ( 656/1549- 698/1886) -43.8 E61-E36 N13-N38
      2 1069 1612 (N29E16 L=119)   375    108        1 (1010/1583-1117/1637)   7.1 E19-E11 N27-N30
      3 1199 1395 (N19E01 L=134)   401    120        2 (1165/1340-1234/1441) -45.0 E03-W02 N16-N21
      4  759 1329 (N03E23 L=112)   242     63        1 ( 743/1298- 775/1366) -76.0 E24-E21 N01-N05
      5 1772 1334 (N36W42 L=176)   384     69        0 (1746/1313-1793/1357)  53.1 W39-W43 N33-N37
      6 1187  386 (S37W34 L=169)   482    134        0 (1129/ 355-1238/ 421)   5.7 W28-W38 S37-S35
      7 1459  335 (S30W55 L=189)   649    155        2 (1423/ 298-1493/ 374)  -5.7 W49-W59 S31-S28
  

    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