Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1653 1552 (N16W60 L=140)   476     61        1 (1622/1533-1682/1570)  -8.1 W57-W64 N14-N16
      2 1587 1449 (N13W48 L=128)  1322    201        1 (1523/1362-1652/1510) -74.7 W45-W50 N06-N17
      3  242 1462 (N53E67 L= 13)   950     94        1 ( 193/1444- 286/1477)   0.0 E86-E56 N51-N53
      4  866 1373 (N29W00 L= 80)  1110    296        3 ( 744/1330- 990/1419) -34.7 E06-W07 N24-N34
      5 1236 1324 (N16W20 L=100)   721    352        0 (1087/1304-1396/1347) -20.2 W11-W30 N12-N20
      6  312 1115 (N28E46 L= 33)   286     52        0 ( 286/1100- 345/1130) -26.6 E49-E43 N26-N29
      7  251  936 (N18E55 L= 24)  1948    395        3 ( 224/ 793- 298/1125)  76.0 E59-E50 N09-N29
      8 1875  997 (S21W60 L=139)   960     76        3 (1847/ 959-1898/1046)  59.0 W55-W63 S24-S18
      9 1366  996 (S05W18 L= 98)   207     45        1 (1347/ 967-1381/1023)  56.3 W17-W19 S06-S03
     10 1675  852 (S23W36 L=116)   548    112        0 (1647/ 807-1704/ 907)  33.7 W33-W39 S24-S20
     11  528  325 (S24E60 L= 20)  1330    125        3 ( 499/ 273- 565/ 397)  25.2 E70-E53 S27-S19
     12  780  285 (S33E42 L= 37)   400     56        0 ( 763/ 255- 798/ 314)  29.7 E45-E38 S35-S31
  

    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