Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  443 1657 (N16E60 L=  9)   719    109        2 ( 397/1631- 496/1680)  20.6 E64-E55 N13-N18
      2  554 1633 (N16E49 L= 20)   358    131        3 ( 525/1607- 606/1677) -11.3 E53-E43 N15-N19
      3 1655 1464 (N38W30 L= 99)   846    160        2 (1632/1384-1682/1547)  77.0 W28-W33 N32-N45
      4 1272 1483 (N26W02 L= 70)   822    301        1 (1162/1371-1348/1545) -19.7 E06-W08 N21-N28
      5  672 1491 (N10E33 L= 35)   287     98        1 ( 655/1442- 693/1530) -68.2 E34-E31 N07-N12
      6  510 1446 (N04E42 L= 26)   294     68        1 ( 487/1425- 541/1467)  56.3 E44-E41 N02-N06
      7  347 1427 (N00E55 L= 14)   955    160        1 ( 271/1393- 414/1458)  20.0 E61-E49 S03-N02
      8  426 1267 (S08E43 L= 26)   350     79        1 ( 400/1236- 452/1298)  56.3 E45-E40 S09-S05
      9  773 1155 (S06E17 L= 52)   236    101        0 ( 755/1108- 786/1205) -80.5 E18-E16 S09-S03
     10 1689 1124 (N18W37 L=106)   360     81        2 (1650/1098-1734/1149)  23.2 W34-W41 N15-N20
     11  257 1054 (S23E53 L= 16)  1911    466        2 (  75/ 993- 430/1136)  14.3 E85-E38 S28-S15
     12 1780 1035 (N16W47 L=116)   350    147        0 (1755/ 963-1795/1098) -45.0 W43-W50 N12-N18
     13  626 1005 (S18E23 L= 46)   858    328        3 ( 459/ 973- 743/1044)  22.8 E34-E15 S21-S13
     14 1714  919 (N07W44 L=113)   338     74        0 (1674/ 902-1750/ 940)  33.7 W41-W47 N05-N09
     15  212  914 (S32E57 L= 12)   408     92        1 ( 181/ 878- 251/ 939)  24.0 E61-E52 S34-S30
     16 1800  776 (N03W57 L=126)   419     79        0 (1785/ 731-1816/ 812) -68.2 W55-W59 N00-N05
     17 1273  667 (S19W25 L= 94)   213     93        0 (1219/ 656-1318/ 681)  18.4 W21-W27 S20-S17
     18 1620  450 (S19W60 L=129)   691    123        3 (1578/ 407-1656/ 483)  33.7 W56-W63 S21-S16
  

    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