Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  369 1624 (N13E67 L=169)  1835    205        5 ( 316/1578- 415/1679) -81.9 E76-E60 N09-N16
      2  456 1585 (N12E55 L=182)  1922    316        4 ( 402/1445- 521/1676) -75.1 E61-E49 N02-N17
      3  562 1631 (N16E49 L=188)   639    134        1 ( 511/1588- 605/1675)  63.4 E51-E46 N13-N19
      4  686 1611 (N18E39 L=198)   297     96        1 ( 676/1562- 705/1658) -54.5 E41-E36 N14-N21
      5 1231 1550 (N29E03 L=233)   333     94        1 (1203/1503-1249/1594) -45.0 E06-W00 N26-N31
      6  512 1548 (N10E48 L=188)   324     64        0 ( 492/1525- 534/1581) -31.0 E50-E45 N09-N12
      7  626 1531 (N12E39 L=198)   425     57        1 ( 599/1497- 651/1560) -21.8 E41-E36 N09-N13
      8 1651 1464 (N39W30 L=267)   654    200        0 (1583/1399-1700/1508) -14.0 W23-W35 N36-N41
      9 1613 1329 (N28W28 L=264)   606    146        0 (1589/1276-1655/1394) -45.0 W25-W32 N26-N32
     10  334 1174 (S15E49 L=188)   279     41        0 ( 314/1156- 358/1194)  36.9 E50-E46 S16-S13
     11  240 1098 (S21E56 L=180)   493    111        1 ( 219/1049- 262/1153)  90.0 E58-E54 S24-S17
     12 1886  959 (N17W62 L=299)   615     82        0 (1848/ 931-1916/ 994)  36.9 W57-W66 N13-N19
     13  487  914 (S27E31 L=206)   311     79        1 ( 445/ 892- 521/ 940)   0.0 E34-E27 S27-S25
     14 1869  801 (N08W66 L=302)   562     48        0 (1857/ 771-1883/ 833) -15.9 W62-W69 N06-N09
     15 1669  600 (S10W55 L=292)   663    143        1 (1613/ 552-1715/ 660)  82.9 W52-W56 S13-S05
     16 1103  589 (S29W18 L=255)   392     96        0 (1054/ 572-1157/ 609)  18.4 W15-W22 S30-S27
  

    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