Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  640 1645 (N35E31 L=308)   402     72        1 ( 599/1623- 679/1661)   0.0 E34-E27 N33-N36
      2  942 1390 (N18E06 L=333)   425     99        2 ( 888/1371- 986/1409)   0.0 E09-E03 N16-N19
      3 1401 1238 (N10W23 L=  1)   533    162        0 (1353/1181-1454/1319) -53.1 W19-W26 N06-N14
      4  932 1152 (N03E06 L=333)   171     51        0 ( 915/1120- 956/1176)  56.3 E06-E04 N01-N04
      5 1604 1057 (S00W37 L= 15)   850    245        0 (1499/1001-1703/1097)  11.3 W29-W44 S03-N02
      6 1828  890 (S09W57 L= 36)  1673    305        2 (1773/ 765-1867/1036)  86.4 W54-W61 S16-N00
      7  522  872 (S13E32 L=307)   267     86        1 ( 506/ 826- 534/ 919)  80.5 E32-E30 S16-S10
      8 1580  857 (S12W36 L= 15)   630    192        5 (1518/ 796-1620/ 898)  -7.1 W31-W39 S15-S09
      9 1869  816 (S12W64 L= 42)  1437    143        4 (1831/ 763-1907/ 883)  54.5 W60-W68 S15-S08
     10  250  836 (S15E55 L=283)   297     81        0 ( 224/ 817- 283/ 869) -36.9 E57-E51 S16-S12
     11  504  619 (S29E37 L=302)  1938    515        4 ( 391/ 524- 667/ 781) -19.3 E45-E24 S35-S18
     12 1569  751 (S18W37 L= 15)   227     54        2 (1541/ 734-1597/ 770)   0.0 W34-W39 S19-S17
     13  313  651 (S26E54 L=285)   551    122        1 ( 291/ 592- 338/ 717)  76.0 E55-E51 S30-S22
     14 1133  654 (S26W08 L=346)   814    196        1 (1038/ 626-1224/ 705)  22.6 W01-W13 S27-S22
     15 1620  550 (S31W47 L= 25)  1655    304        3 (1506/ 490-1738/ 652)  32.0 W37-W55 S34-S23
     16 1429  467 (S37W33 L= 12)   291     49        1 (1410/ 446-1449/ 493)  71.6 W31-W35 S38-S35
     17 1153  424 (S42W12 L=350)  1741    397        1 ( 963/ 398-1336/ 469)  10.1 E03-W25 S43-S37
  

    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