Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  915 1698 (N30E28 L=326)   333     87        1 ( 864/1682- 963/1710)  38.7 E31-E25 N27-N31
      2 1665 1453 (N39W31 L= 25)  3438    514        4 (1520/1276-1761/1653) -42.3 W15-W40 N28-N49
      3 1465 1630 (N44W11 L=  5)   657    107        0 (1422/1601-1502/1649)  29.7 W07-W14 N41-N47
      4  385 1488 (N04E55 L=299)   536    109        0 ( 348/1449- 409/1513)  -8.1 E59-E52 N01-N05
      5 1033 1268 (N07E06 L=348)   257     83        0 (1018/1224-1054/1309) -51.3 E07-E03 N04-N09
      6  167 1277 (S12E70 L=285)   557    186        0 ( 134/1260- 208/1295)  11.3 E74-E64 S12-S10
      7 1750 1246 (N28W40 L= 35)   766    112        0 (1716/1207-1771/1293)  63.4 W37-W42 N25-N31
      8 1538 1171 (N16W26 L= 20)   420    119        0 (1490/1110-1565/1234) -35.5 W21-W28 N13-N18
      9 1797 1105 (N21W47 L= 42)   460     99        0 (1764/1081-1831/1134)   8.1 W43-W50 N20-N23
     10 1459 1019 (N05W24 L= 19)   917    210        2 (1392/ 968-1555/1098)  -4.4 W18-W31 N03-N07
     11 1868  975 (N17W58 L= 52)   916    143        2 (1846/ 908-1884/1038) -50.2 W54-W61 N13-N20
     12 1704  945 (N09W43 L= 38)  1233    258        2 (1612/ 914-1790/ 970)  26.6 W37-W51 N04-N12
     13 1714  763 (S00W51 L= 45)  1024    176        0 (1642/ 700-1771/ 842)  68.2 W47-W53 S05-N05
     14  724  757 (S30E11 L=344)   234     59        0 ( 691/ 734- 747/ 787)  53.1 E12-E09 S31-S27
     15 1125  644 (S25W17 L= 12)  1179    305        2 ( 990/ 533-1257/ 721)  -5.0 W06-W29 S27-S21
     16 1515  599 (S15W43 L= 37)   446    116        0 (1459/ 563-1565/ 635)  54.5 W40-W45 S18-S11
  

    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