Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  431 1597 (N20E59 L=256)   473     70        0 ( 410/1559- 455/1628) -78.7 E60-E57 N17-N22
      2  658 1575 (N23E38 L=277)   276     70        1 ( 639/1552- 686/1604) -11.3 E40-E35 N21-N25
      3  979 1527 (N27E15 L=300)   980    296        3 ( 849/1452-1104/1583)  -5.4 E25-E04 N25-N29
      4 1228 1423 (N26W04 L=319)   739    351        0 (1131/1311-1319/1556) -36.9 E05-W11 N21-N33
      5  539 1419 (N11E40 L=275)   303     50        1 ( 520/1391- 556/1445) -56.3 E41-E38 N09-N12
      6 1751 1230 (N27W47 L=  3)  1493    237        2 (1717/1129-1796/1330)  70.3 W43-W54 N20-N34
      7  276 1138 (S10E54 L=262)   494     59        0 ( 244/1120- 306/1155)  18.4 E57-E51 S11-S08
      8 1579 1048 (N11W34 L=349)   284     44        2 (1556/1021-1603/1068)  45.0 W32-W35 N09-N13
      9 1874  957 (N13W65 L= 20)  1812    110        5 (1843/ 912-1920/1007)  80.5 W59-W73 N10-N16
     10  327  930 (S22E46 L=269)   727    157        4 ( 287/ 881- 361/ 970) -78.7 E50-E42 S24-S19
     11 1797  920 (N09W56 L= 11)   423     72        1 (1784/ 880-1813/ 952) -63.4 W53-W57 N06-N11
     12  404  874 (S24E39 L=277)   698    167        6 ( 376/ 799- 432/ 947) -71.6 E41-E36 S28-S19
     13 1236  631 (S20W22 L=338)  1187    272        3 (1183/ 517-1337/ 720)  82.9 W20-W25 S28-S12
     14  327  591 (S45E50 L=265)  1085    152        3 ( 306/ 508- 347/ 651)  65.6 E56-E46 S51-S40
     15  723  598 (S35E10 L=305)   350    124        2 ( 668/ 567- 786/ 631)  32.0 E14-E06 S37-S30
     16 1317  497 (S27W33 L=348)   265    118        1 (1254/ 481-1367/ 509)  14.0 W28-W36 S27-S25
     17 1517  408 (S27W54 L= 10)   434     81        0 (1501/ 370-1538/ 442) -45.0 W51-W56 S30-S25
  

    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