Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  607 1740 (N31E49 L= 65)  1352    196        3 ( 527/1677- 679/1794) -13.5 E64-E38 N26-N34
      2  968 1690 (N33E19 L= 94)   642    110        1 ( 931/1633-1019/1731) -29.1 E22-E13 N30-N36
      3 1076 1497 (N22E07 L=107)   528    142        2 (1032/1448-1130/1560) -26.6 E10-E02 N19-N25
      4  905 1332 (N09E13 L=101)  1277    438        2 ( 876/1169- 935/1542) -62.4 E18-E07 N00-N21
      5 1640 1348 (N26W34 L=148)   655    133        0 (1609/1298-1680/1407) -45.0 W31-W37 N23-N29
      6 1236 1364 (N17W06 L=120)   393     73        3 (1207/1330-1274/1396)  59.0 W04-W07 N14-N19
      7 1860 1147 (N20W57 L=171)   531     53        0 (1829/1121-1886/1171)  -8.1 W53-W60 N18-N21
      8 1248  847 (S12W17 L=130)   376    160        1 (1177/ 814-1320/ 880)  32.0 W12-W20 S15-S09
      9 1062  800 (S19W07 L=121)   253     53        1 (1042/ 766-1087/ 825) -56.3 W05-W08 S20-S17
     10  775  783 (S25E11 L=103)   357    100        0 ( 729/ 747- 820/ 816) -18.4 E14-E07 S26-S23
     11  711  650 (S35E13 L=101)   317     87        1 ( 680/ 630- 755/ 682)  -9.5 E15-E09 S36-S33
     12 1715  643 (S12W56 L=170)   439     56        1 (1691/ 624-1742/ 667)   0.0 W53-W59 S13-S11
     13 1473  460 (S29W46 L=160)  1101    293        0 (1333/ 394-1564/ 496)   4.8 W34-W59 S31-S26
     14  831  455 (S46E00 L=114)   398     59        1 ( 801/ 431- 863/ 485) -31.0 E02-W03 S47-S43
     15 1317  238 (S47W54 L=168)   483     74        1 (1278/ 218-1357/ 256)  11.3 W49-W59 S48-S46
     16  800  215 (S68W07 L=121)   826    131        2 ( 727/ 187- 859/ 239)   2.7 E05-W17 S69-S65
  

    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