Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1205 1396 (N27W13 L=196)   403    113        2 (1178/1330-1219/1442) -90.0 W11-W14 N22-N30
      2  393 1328 (N23E45 L=137)   384     84        1 ( 377/1282- 411/1365) -80.5 E46-E44 N19-N25
      3 1604 1277 (N18W41 L=224)   899    156        3 (1553/1242-1664/1308) -26.6 W36-W45 N15-N19
      4 1333 1165 (N12W20 L=203)   197     99        1 (1277/1153-1374/1175)  -9.5 W16-W22 N10-N12
      5 1770  837 (S11W53 L=236)  1645    371       10 (1720/ 755-1826/ 941)  70.0 W48-W59 S15-S04
      6 1676  777 (S14W45 L=228)   400     84        3 (1650/ 738-1698/ 815)  59.0 W42-W46 S16-S11
      7 1745  667 (S22W53 L=236)  2097    547        3 (1594/ 588-1880/ 761)  18.4 W40-W69 S26-S16
      8 1418  696 (S18W25 L=208)   371    144        1 (1374/ 636-1453/ 739) -54.5 W22-W28 S21-S14
      9  510  676 (S18E36 L=147)   257     53        1 ( 495/ 649- 529/ 699)  71.6 E37-E34 S19-S16
     10  863  466 (S33E13 L=170)   271     84        0 ( 819/ 447- 909/ 479)   8.1 E16-E09 S33-S31
     11 1630  422 (S39W51 L=234)   799    152        4 (1575/ 396-1674/ 453)  14.0 W45-W54 S41-S37
  

    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