Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  202 1158 (S11E61 L=248)   382     79        1 ( 173/1134- 230/1183)  -8.1 E64-E57 S12-S09
      2  329 1125 (S12E47 L=261)   715    181        2 ( 259/1078- 397/1163)  31.0 E52-E42 S15-S08
      3  723 1037 (S11E18 L=291)   447    115        1 ( 694/ 997- 758/1094) -51.3 E19-E15 S13-S08
      4 1867 1054 (N14W58 L=  6)   500    176        1 (1822/1030-1904/1072)  24.0 W53-W62 N11-N15
      5  617  983 (S16E24 L=285)   585    179        1 ( 538/ 941- 706/1042)  -9.5 E29-E17 S17-S13
      6  377  909 (S24E41 L=268)  1002    337        1 ( 272/ 834- 499/1031)  39.5 E49-E32 S29-S15
      7  161  950 (S24E62 L=247)  1592    283        2 ( 124/ 874- 194/1024)  68.2 E65-E58 S28-S18
      8 1855  928 (N06W59 L=  8)  1240    199        4 (1823/ 886-1884/ 980)  80.5 W56-W63 N03-N09
      9 1145  880 (S13W10 L=318)   180     53        0 (1111/ 867-1178/ 894)  33.7 W08-W11 S14-S11
     10 1039  658 (S28W09 L=317)   258     77        0 ( 998/ 629-1075/ 688)  51.3 W06-W10 S30-S25
     11 1598  653 (S15W45 L=354)   332     66        0 (1566/ 632-1628/ 676)  53.1 W43-W47 S17-S13
     12  527  626 (S40E28 L=281)   389    106        1 ( 491/ 597- 574/ 660)  35.5 E31-E23 S42-S36
     13  322  579 (S46E49 L=260)   513     52        0 ( 307/ 546- 337/ 613) -90.0 E50-E47 S48-S43
     14 1644  565 (S19W53 L=  1)   328     88        0 (1604/ 551-1677/ 582)  31.0 W49-W55 S20-S17
     15 1608  488 (S24W54 L=  2)   462    128        0 (1541/ 461-1659/ 513)  40.6 W49-W57 S27-S21
     16 1562  338 (S34W59 L=  8)  1737    328        5 (1467/ 306-1641/ 370)  21.8 W51-W66 S36-S30
  

    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