Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1461 1730 (N57W36 L=167)  1325    199        3 (1385/1683-1554/1753)  14.0 W23-W55 N51-N60
      2  921 1645 (N40E18 L=114)   323     69        1 ( 882/1623- 955/1661)   9.5 E21-E15 N38-N41
      3 1439 1626 (N47W27 L=159)   326     48        1 (1407/1610-1464/1647)  18.4 W24-W30 N45-N48
      4  677 1446 (N23E29 L=102)   571    144        2 ( 646/1388- 711/1516) -71.6 E32-E26 N18-N27
      5 1221 1461 (N31W08 L=139)   300     85        2 (1199/1420-1235/1506) -68.2 W06-W09 N28-N33
      6  482 1428 (N19E44 L= 88)   589     85        4 ( 452/1392- 520/1464)  59.0 E47-E40 N16-N21
      7 1748 1413 (N34W56 L=188)   819    104        1 (1702/1386-1795/1435) -10.3 W50-W62 N32-N35
      8  469  814 (S19E35 L= 97)   788    157        2 ( 426/ 746- 525/ 854)  41.2 E38-E30 S23-S15
      9 1426  686 (S15W31 L=162)   252     94        1 (1393/ 649-1457/ 720)  53.1 W29-W32 S17-S12
     10 1710  649 (S14W56 L=187)   445     39        2 (1690/ 622-1730/ 676) -63.4 W53-W57 S16-S12
     11 1381  385 (S36W39 L=170)   397    100        0 (1343/ 359-1426/ 424)  63.4 W36-W40 S38-S32
  

    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