Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  649 1664 (N42E39 L=168)  1061    153        2 ( 616/1583- 679/1724) -38.2 E46-E32 N35-N46
      2 1079 1483 (N32W00 L=208)  1559    445        1 ( 950/1334-1232/1571) -33.0 E09-W11 N22-N37
      3  405 1390 (N20E49 L=159)   616    146        1 ( 386/1334- 436/1470) -48.4 E53-E44 N16-N25
      4 1597 1330 (N25W39 L=247)   272     39        0 (1580/1309-1619/1348) -18.4 W37-W41 N23-N25
      5 1774 1168 (N15W54 L=262)  1097    193        1 (1738/1091-1805/1246) -90.0 W51-W57 N10-N19
      6 1573 1161 (N14W36 L=243)   238     53        1 (1552/1137-1596/1187) -33.7 W34-W37 N12-N15
      7 1843 1077 (N10W62 L=269)   626     92        1 (1815/1039-1868/1116)  29.7 W58-W65 N07-N12
      8 1754  858 (S04W53 L=260)   520     11        0 (1730/ 834-1777/ 881) -36.9 W50-W55 S05-S02
      9  972  754 (S15E01 L=206)   242     56        2 ( 954/ 719- 990/ 781)  90.0 E02-W00 S17-S13
     10 1856  715 (S13W72 L=279)  1844    164        5 (1826/ 648-1886/ 768) -82.9 W66-W78 S17-S09
     11 1176  671 (S19W12 L=219)   255     80        0 (1150/ 633-1200/ 711)  68.2 W10-W13 S21-S16
     12  358  672 (S25E48 L=159)   339     57        1 ( 336/ 648- 379/ 701)  45.0 E50-E45 S27-S23
     13 1755  647 (S18W60 L=267)   708     42        1 (1726/ 624-1783/ 670) -12.5 W55-W64 S19-S16
     14  282  585 (S33E62 L=145)   535     55        0 ( 249/ 574- 315/ 597)   0.0 E68-E57 S33-S31
     15 1226  577 (S25W17 L=224)   284     77        0 (1186/ 563-1264/ 594)  21.8 W14-W19 S25-S23
     16  904  494 (S33E05 L=202)   281     67        0 ( 869/ 479- 941/ 510)  31.0 E07-E02 S34-S31
     17 1089  425 (S37W10 L=217)   378    116        0 (1035/ 406-1148/ 443)  18.4 W05-W14 S39-S35
  

    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