Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1606 1716 (N28W74 L= 44)   903     91        0 (1579/1687-1641/1748)  59.0 W69-W78 N25-N30
      2 1528 1665 (N29W59 L= 28)   604     25        0 (1507/1641-1552/1694)  56.3 W55-W62 N27-N30
      3  573 1639 (N58E21 L=308)   320     53        0 ( 550/1617- 601/1659)  18.4 E24-E17 N56-N59
      4 1076 1579 (N38W20 L=350)   477    101        1 (1031/1540-1110/1603) -63.4 W17-W23 N34-N40
      5 1651 1551 (N19W60 L= 30)  2184    342        3 (1491/1512-1783/1586) -15.7 W45-W77 N13-N23
      6  297 1386 (N45E52 L=277)   832     93        0 ( 255/1366- 348/1405) -14.9 E59-E44 N42-N46
      7  497 1321 (N36E29 L=301)   298     59        1 ( 470/1296- 526/1342)  11.3 E31-E26 N34-N37
      8 1387  845 (S12W17 L=346)  1118    295        1 (1339/ 723-1421/ 990)  52.4 W11-W21 S17-S03
      9  733  714 (S04E24 L=306)   452    144        1 ( 710/ 643- 762/ 778)  45.0 E27-E20 S07-N00
     10 1062  387 (S32E15 L=315)  2045    344        6 ( 911/ 283-1154/ 496)  23.5 E29-E06 S36-S26
     11  562  418 (S19E47 L=283)   549     65        1 ( 529/ 394- 589/ 443)   9.5 E50-E43 S20-S17
     12  777  397 (S25E33 L=297)   365     87        0 ( 744/ 364- 804/ 437)  29.7 E36-E29 S26-S22
  

    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