Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  742 1825 (N62E42 L= 12)   811    116        1 ( 696/1807- 794/1844)   8.1 E51-E31 N59-N63
      2 1053 1709 (N51W02 L= 55)   320     68        1 (1030/1676-1072/1729)  53.1 W00-W03 N48-N52
      3  551 1615 (N42E44 L=  9)   367     69        0 ( 540/1573- 569/1649) -50.2 E46-E41 N38-N44
      4  952 1368 (N25E05 L= 48)   413    177        2 ( 903/1305-1015/1424) -41.2 E08-W00 N21-N28
      5 1705 1396 (N26W54 L=107)   458     97        1 (1673/1372-1756/1411)  -6.3 W50-W60 N24-N27
      6 1681  850 (S08W46 L= 99)   794    216        2 (1583/ 827-1767/ 877)  -3.6 W37-W53 S09-S06
      7 1812  838 (S09W60 L=113)   535     81        0 (1789/ 810-1831/ 873)  56.3 W56-W61 S11-S07
      8 1199  732 (S15W11 L= 65)   887    297        3 (1114/ 680-1278/ 810) -16.7 W05-W16 S18-S09
      9  596  690 (S18E29 L= 24)   960    176        2 ( 529/ 627- 669/ 729)  28.6 E34-E23 S22-S15
     10 1123  628 (S22W07 L= 60)   240     78        1 (1085/ 605-1162/ 651)  21.8 W04-W09 S23-S20
     11 1245  575 (S25W16 L= 69)  1094    439        2 (1067/ 533-1438/ 623)  10.9 W03-W29 S28-S22
     12 1561  544 (S28W41 L= 95)   432     88        0 (1521/ 516-1596/ 578)  38.7 W38-W43 S29-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