Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1061 1715 (N36E15 L= 89)   761    145        2 ( 994/1670-1115/1748) -16.7 E20-E10 N34-N38
      2  419 1562 (N14E53 L= 51)   339     59        1 ( 388/1545- 451/1577)  21.8 E56-E50 N12-N15
      3 1324 1511 (N28W08 L=112)   726    247        0 (1221/1466-1424/1542)  30.3 W02-W14 N23-N31
      4 1587 1499 (N35W28 L=132)  1071    200        2 (1502/1441-1665/1534)   4.1 W20-W34 N32-N37
      5 1117 1250 (N08W01 L=105)   312     77        0 (1079/1229-1162/1269)  36.9 E01-W03 N05-N09
      6  216 1222 (S09E60 L= 44)   500     86        1 ( 189/1172- 237/1262)  71.6 E62-E58 S12-S06
      7  345 1132 (S13E46 L= 58)   485    115        1 ( 289/1104- 401/1158)  -5.7 E50-E40 S13-S11
      8  616 1095 (S11E25 L= 79)   248     96        0 ( 596/1048- 632/1144) -80.5 E26-E24 S13-S07
      9 1845 1086 (N17W54 L=158)   896     95        0 (1800/1042-1886/1134) -11.3 W48-W59 N14-N18
     10  923  964 (S13E04 L=100)   937    491        0 ( 714/ 881-1116/1041)  35.0 E15-W05 S21-S05
     11 1835  748 (S02W63 L=167)  2530    514        4 (1742/ 548-1884/ 969) -78.7 W58-W72 S15-N10
     12  326  948 (S24E46 L= 59)   331     59        1 ( 290/ 930- 363/ 966)  18.4 E48-E42 S25-S22
     13 1419  865 (S07W26 L=130)   210     56        0 (1386/ 849-1450/ 880)   0.0 W24-W28 S07-S06
     14 1158  825 (S15W12 L=116)   245     91        0 (1134/ 783-1183/ 858)  78.7 W10-W13 S18-S13
     15 1243  808 (S14W17 L=121)   232    116        1 (1210/ 773-1279/ 840)  53.1 W15-W19 S16-S11
     16 1232  685 (S21W20 L=124)   598    241        1 (1125/ 652-1345/ 737)  24.8 W13-W26 S24-S16
     17 1176  584 (S29W20 L=124)   296     85        0 (1140/ 565-1216/ 603)  45.0 W17-W22 S30-S26
  

    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