Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  431 1637 (N14E60 L=281)   488     83        0 ( 416/1604- 445/1671) -53.1 E61-E57 N12-N16
      2 1118 1629 (N31E13 L=328)   251     50        1 (1083/1616-1147/1645)  36.9 E15-E11 N29-N32
      3 1594 1572 (N46W23 L=  4)   752    142        0 (1565/1528-1617/1618)  73.3 W20-W25 N40-N50
      4  528 1446 (N05E42 L=299)   320    111        0 ( 482/1413- 563/1467)   0.0 E45-E38 N02-N05
      5 1769 1423 (N42W42 L= 23)  1308    161        1 (1716/1384-1821/1451)  10.3 W35-W48 N37-N44
      6 1848 1263 (N34W50 L= 31)  2665    345        6 (1778/1127-1879/1358) -58.4 W42-W54 N25-N40
      7  202 1240 (S13E63 L=278)  2163    475        3 ( 109/1202- 291/1269)   0.0 E77-E52 S14-S11
      8 1234 1159 (N06W08 L=349)   253    132        0 (1193/1130-1260/1217) -53.1 W06-W09 N03-N09
      9 1700 1074 (N16W39 L= 20)   335    102        0 (1667/1021-1720/1117) -38.7 W36-W41 N13-N17
     10  483 1031 (S20E33 L=308)   229     70        0 ( 448/1009- 513/1049)   0.0 E35-E30 S21-S18
     11 1621  925 (N05W38 L= 19)   641    158        0 (1562/ 875-1679/ 974)  -5.7 W32-W43 N03-N06
     12 1836  864 (N10W58 L= 39)  1554    310        5 (1763/ 826-1902/ 897)  17.4 W50-W67 N05-N13
     13 1446  688 (S12W34 L= 16)   476    150        1 (1392/ 652-1501/ 752)  60.3 W31-W36 S15-S07
     14 1728  601 (S07W60 L= 41)  2072    394        1 (1588/ 488-1852/ 742)  60.3 W53-W67 S17-N04
     15  831  653 (S33E01 L=340)   489    185        1 ( 747/ 629- 927/ 696)  45.0 E05-W04 S36-S27
     16 1226  567 (S26W26 L=  7)   844    177        2 (1137/ 520-1297/ 632)  -4.4 W18-W31 S27-S24
  

    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