Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1221 1844 (N58W47 L=119)   425     60        0 (1205/1815-1237/1873)  22.6 W41-W53 N55-N60
      2  680 1717 (N59E23 L= 49)   402     84        0 ( 644/1697- 724/1736)  12.5 E26-E16 N56-N60
      3 1629 1634 (N32W65 L=137)   746     98        0 (1600/1591-1652/1682)  13.2 W58-W75 N30-N34
      4 1126 1513 (N35W15 L= 87)   285     67        0 (1088/1498-1166/1527)   0.0 W12-W18 N33-N35
      5  902 1376 (N29E03 L= 68)   487    179        0 ( 842/1302- 937/1427) -66.0 E06-E01 N24-N33
      6  686 1338 (N30E19 L= 53)   320    125        0 ( 651/1305- 739/1389) -54.5 E21-E14 N26-N33
      7  157 1191 (N24E70 L=  1)   792     97        0 ( 142/1152- 177/1242) -45.0 E75-E66 N21-N27
      8 1397 1191 (N11W25 L= 97)   200     28        1 (1379/1170-1413/1216)  56.3 W24-W26 N09-N12
      9  282 1152 (N21E53 L= 19)   457    107        2 ( 266/1104- 301/1203) -71.6 E54-E50 N18-N24
     10 1825  956 (S11W56 L=128)  1885    281        2 (1780/ 870-1871/1024)  45.0 W50-W62 S15-S07
     11  658  872 (N01E24 L= 47)   328     80        1 ( 628/ 841- 689/ 907) -45.0 E26-E22 N00-N03
     12 1714  622 (S31W46 L=118)  2569    483        5 (1540/ 491-1876/ 746)  12.5 W30-W67 S37-S25
     13 1308  390 (S40W11 L= 83)   499     97        0 (1252/ 369-1353/ 416)  -6.3 W06-W15 S41-S38
     14 1142  332 (S42E03 L= 69)   448     67        1 (1112/ 311-1183/ 352)  11.3 E05-W00 S43-S39
     15  963  246 (S47E21 L= 51)  1228    144        2 ( 898/ 214-1035/ 282)   4.1 E28-E13 S48-S43
     16  809  191 (S50E40 L= 32)   812     88        0 ( 785/ 148- 835/ 232)  29.1 E51-E33 S55-S45
  

    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