Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  462 1664 (N58E53 L= 29)   838     89        0 ( 427/1636- 509/1686)  18.4 E60-E45 N56-N60
      2 1619 1436 (N20W50 L=133)   424    146        0 (1564/1398-1668/1460) -26.6 W45-W53 N16-N21
      3  649 1411 (N36E21 L= 61)   390    107        0 ( 590/1394- 705/1427)  -6.3 E25-E16 N34-N37
      4  292 1060 (N17E52 L= 30)   531     87        0 ( 274/1027- 314/1105) -59.0 E53-E49 N15-N20
      5 1885 1010 (S11W65 L=148)   757     74        1 (1851/ 986-1923/1036)   0.0 W60-W72 S12-S10
      6 1381  676 (S21W17 L=100)   315     96        1 (1343/ 633-1416/ 721)  33.7 W14-W20 S23-S18
      7 1070  575 (S22E04 L= 78)   941    271        1 ( 983/ 497-1158/ 681)  37.6 E11-W02 S26-S16
      8 1539  572 (S32W29 L=112)   393    133        0 (1489/ 522-1581/ 626)  40.6 W25-W32 S34-S28
      9  585  568 (S16E37 L= 46)   403     95        0 ( 561/ 522- 622/ 604)  33.7 E39-E33 S18-S13
     10 1577  512 (S37W34 L=116)   408    101        0 (1540/ 492-1618/ 549)  15.9 W30-W37 S38-S35
     11 1108  465 (S31E04 L= 78)   927    348        0 ( 971/ 398-1249/ 541)  13.4 E14-W07 S33-S27
     12  695  428 (S27E34 L= 48)   498     88        0 ( 675/ 391- 718/ 469)  80.5 E35-E31 S29-S23
  

    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