Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  769 1822 (N62E37 L= 13)   713    111        0 ( 724/1801- 818/1840)  -6.0 E46-E27 N59-N63
      2 1032 1354 (N24W00 L= 50)   850    256        1 ( 953/1278-1144/1424) -26.6 E05-W07 N19-N28
      3 1736 1400 (N26W58 L=107)   553     99        0 (1701/1376-1787/1416)  -5.7 W53-W65 N24-N27
      4 1717  853 (S08W49 L= 98)  1037    197        1 (1630/ 828-1788/ 884)  -3.8 W41-W56 S09-S05
      5 1206  754 (S13W12 L= 61)   357    134        0 (1172/ 723-1255/ 806) -36.9 W09-W14 S15-S09
      6 1309  740 (S14W19 L= 68)   401    101        1 (1274/ 701-1340/ 776)  45.0 W16-W20 S16-S11
      7  646  686 (S18E25 L= 24)  1081    192        2 ( 576/ 624- 728/ 729)  30.3 E31-E19 S22-S15
      8  344  535 (S30E57 L=352)  1331    178        3 ( 323/ 459- 373/ 647)  51.8 E65-E51 S36-S22
      9  258  593 (S27E67 L=342)  1500    102        2 ( 228/ 539- 278/ 633)  37.9 E72-E61 S30-S23
     10 1252  577 (S25W16 L= 65)  1481    536        3 (1104/ 533-1473/ 628)  13.5 W05-W31 S28-S21
     11 1600  548 (S28W45 L= 94)   498    123        0 (1558/ 507-1649/ 598)  45.0 W41-W48 S30-S24
  

    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