Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1257 1763 (N46E09 L=102)   305     39        1 (1232/1745-1284/1778)  53.1 E10-E07 N44-N48
      2 1187 1694 (N38E11 L=101)  1124    244        2 (1125/1617-1302/1738)  65.0 E13-E04 N30-N45
      3  786 1650 (N23E33 L= 78)   260     29        0 ( 764/1628- 807/1669) -26.6 E35-E31 N21-N24
      4 1634 1454 (N37W28 L=140)   404     53        0 (1613/1420-1652/1477) -56.3 W26-W30 N34-N38
      5 1581 1270 (N22W27 L=138)   322     63        1 (1559/1232-1600/1306)  80.5 W25-W27 N19-N25
      6 1818 1071 (N19W50 L=162)  2281    481        8 (1769/ 901-1870/1273) -55.6 W42-W56 N10-N30
      7 1184 1231 (N08W03 L=115)   450    111        1 (1135/1200-1246/1257)  39.8 W00-W06 N05-N11
      8  122 1142 (S19E73 L= 39)  1248     84        1 (  77/1122- 159/1163)   5.7 E87-E67 S20-S17
      9  748 1040 (S13E16 L= 96)   200     27        0 ( 729/1019- 767/1062) -18.4 E17-E14 S13-S11
     10  669  958 (S19E19 L= 93)   265    117        0 ( 616/ 943- 732/ 969)  14.0 E22-E14 S21-S18
     11 1860  862 (N10W61 L=172)   841     93        3 (1821/ 825-1899/ 891)  24.0 W56-W66 N07-N12
     12 1093  807 (S17W10 L=121)   404     62        2 (1064/ 778-1121/ 840)  76.0 W07-W11 S18-S14
     13 1527  599 (S15W43 L=155)   382    143        1 (1463/ 574-1579/ 637)  54.5 W40-W45 S18-S11
     14 1640  489 (S17W59 L=171)  1212    150        2 (1587/ 442-1688/ 546)  71.6 W55-W66 S19-S13
     15 1681  405 (S19W72 L=183)  1680    159        3 (1616/ 360-1761/ 449)  31.0 W66-W84 S23-S14
     16 1506  316 (S30W63 L=175)   780     34        1 (1482/ 287-1534/ 349) -12.5 W58-W67 S31-S28
  

    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