Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  478 1658 (N37E58 L=343)  1101    109        0 ( 423/1632- 540/1685)  -2.9 E68-E48 N35-N39
      2  587 1546 (N31E39 L=  1)   716    109        1 ( 559/1489- 633/1593) -49.4 E43-E35 N26-N33
      3  922 1418 (N26E11 L= 30)  1303    341        0 ( 827/1343- 973/1540) -72.9 E16-E06 N20-N33
      4 1747 1375 (N30W56 L= 96)  3558    517        1 (1553/1205-1876/1517) -33.2 W37-W68 N19-N38
      5  350 1373 (N17E53 L=347)   374     66        0 ( 318/1349- 379/1391)   0.0 E56-E50 N15-N18
      6 1263 1311 (N21W13 L= 54)   458    216        0 (1187/1268-1329/1357)  33.7 W08-W17 N18-N24
      7 1784 1196 (N18W55 L= 96)  1422    307        1 (1721/1060-1822/1314) -66.4 W50-W59 N09-N25
      8 1905 1202 (N19W74 L=115)  1377     78        3 (1884/1159-1926/1242)  68.2 W70-W79 N16-N21
      9 1691 1157 (N15W46 L= 86)   577    159        0 (1651/1115-1739/1229) -40.6 W42-W50 N12-N19
     10 1397 1125 (N11W23 L= 63)   233     72        0 (1364/1105-1425/1149)  36.9 W20-W24 N09-N12
     11  808  728 (S19E11 L= 29)   700    225        0 ( 777/ 638- 874/ 786)  73.3 E13-E07 S24-S14
     12 1741  640 (S17W59 L= 99)  1065    157        0 (1688/ 598-1794/ 701)  41.2 W54-W63 S19-S12
  

    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