Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1247 1461 (N28W25 L=183)   597    113        2 (1197/1423-1301/1493)  14.0 W21-W29 N25-N29
      2  444 1272 (N33E35 L=123)  1534    390        2 ( 356/1147- 560/1396) -49.6 E44-E27 N22-N42
      3  192 1213 (N32E60 L= 99)   562     88        1 ( 171/1192- 221/1231)  11.3 E61-E56 N30-N32
      4  227 1066 (N22E56 L=102)  1564    334        3 ( 179/ 918- 280/1163) -86.2 E61-E51 N13-N28
      5 1657 1099 (S03W42 L=200)   957    161        2 (1600/1067-1720/1137) -40.6 W36-W46 S06-N00
      6  671 1065 (N16E21 L=138)   278    114        0 ( 625/1053- 732/1085) -23.2 E23-E16 N14-N17
      7  830  895 (N03E14 L=144)  1078    238        5 ( 769/ 819- 900/ 980) -74.7 E16-E10 S02-N09
      8 1524  824 (S16W27 L=185)   485     59        0 (1507/ 785-1544/ 868)  63.4 W25-W28 S18-S13
      9  812  576 (S16E22 L=136)   823    173        3 ( 727/ 515- 864/ 617) -53.1 E27-E19 S20-S12
     10  679  283 (S34E45 L=113)  1016    263        1 ( 635/ 236- 733/ 338)  25.0 E52-E37 S37-S30
  

    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