Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1455 1683 (N31W50 L=268)   371     84        2 (1422/1647-1483/1715)  78.7 W47-W51 N28-N33
      2  797 1675 (N48W02 L=220)   259     56        0 ( 762/1657- 822/1688) -26.6 E01-W03 N46-N49
      3 1603 1626 (N23W58 L=277)   890    106        3 (1554/1583-1657/1675) -50.2 W55-W62 N20-N27
      4  601 1629 (N51E17 L=201)   484     71        0 ( 565/1602- 633/1663)  18.4 E20-E13 N49-N53
      5 1376 1588 (N27W38 L=256)  3894    776        8 (1070/1496-1749/1648) -20.4 W18-W70 N15-N35
      6 1684 1335 (N05W49 L=268)  1045    172        3 (1599/1287-1742/1370) -26.6 W43-W54 N02-N08
      7  193 1256 (N33E60 L=159)   572     74        1 ( 159/1220- 224/1285) -24.4 E66-E54 N30-N35
      8  618 1174 (N20E21 L=197)  1429    448        3 ( 425/1108- 781/1237) -31.3 E34-E11 N12-N27
      9  216 1128 (N25E55 L=163)   367     56        0 ( 196/1101- 237/1157)  45.0 E57-E52 N23-N26
     10  141 1108 (N25E66 L=153)   662     85        1 ( 111/1083- 182/1135) -14.0 E71-E59 N23-N27
     11  331  942 (N12E46 L=173)   473    116        3 ( 305/ 884- 347/ 998) -81.9 E47-E43 N08-N15
     12 1482  940 (S12W25 L=243)   266     58        1 (1461/ 915-1510/ 973) -63.4 W23-W26 S13-S09
     13 1443  817 (S18W20 L=239)   244     63        0 (1406/ 796-1478/ 832) -36.9 W17-W22 S19-S16
     14  889  722 (S11E14 L=204)   402    127        0 ( 869/ 669- 929/ 774)  56.3 E15-E10 S14-S08
     15 1886  746 (S35W63 L=281)   997    100        0 (1834/ 724-1939/ 765)  -7.4 W54-W77 S36-S33
     16 1500  643 (S30W22 L=240)   910    236        0 (1392/ 587-1567/ 716)  12.1 W13-W27 S32-S26
     17  866  443 (S28E24 L=194)  1864    458        2 ( 683/ 365-1072/ 554)   5.4 E40-E08 S28-S24
     18 1240  341 (S45E03 L=215)  1836    367        2 (1073/ 290-1376/ 443)   6.6 E17-W09 S47-S40
  

    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