Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1308 1700 (N35W45 L=148)  1374    250        1 (1214/1667-1427/1722) -32.7 W36-W52 N29-N38
      2  393 1412 (N46E37 L= 67)  1649    198        3 ( 303/1348- 484/1464) -22.6 E51-E27 N39-N51
      3 1560 1164 (S02W35 L=138)   512    110        0 (1540/1107-1582/1217)  80.5 W33-W36 S05-N01
      4 1832 1051 (S17W54 L=158)  1076    218        1 (1760/ 973-1894/1128)  11.9 W45-W64 S18-S14
      5  159 1068 (N28E66 L= 38)  1709    166        2 (  92/1008- 223/1108) -14.6 E81-E56 N23-N31
      6  142  845 (N15E72 L= 32)   807     92        2 ( 104/ 812- 178/ 877) -26.6 E80-E67 N12-N17
      7 1610  681 (S31W28 L=132)  1402    313        0 (1567/ 559-1640/ 817)  73.6 W22-W30 S38-S21
      8  441  751 (N05E43 L= 61)   278     52        0 ( 409/ 732- 471/ 771) -26.6 E45-E40 N03-N06
      9  253  734 (N07E61 L= 43)   559     67        1 ( 227/ 713- 281/ 760)   0.0 E63-E56 N05-N09
     10  766  437 (S21E33 L= 70)  1027    178        1 ( 702/ 372- 825/ 506)  21.0 E40-E26 S24-S17
     11 1368  416 (S41W04 L=107)   378     79        0 (1325/ 397-1414/ 431) -15.9 W00-W07 S42-S40
  

    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