Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1595 1528 (N26W46 L=  2)   741    132        1 (1531/1490-1654/1557) -21.8 W40-W50 N23-N28
      2 1808  976 (S09W55 L= 11)   874    214        1 (1757/ 895-1863/1065)  42.3 W49-W61 S13-S03
      3 1457 1033 (S04W27 L=343)   241     55        0 (1431/1011-1491/1051)  21.8 W24-W29 S05-S02
      4  925  804 (S14E08 L=309)   451    258        1 ( 842/ 725- 998/ 901) -56.3 E11-E03 S19-S07
      5  465  801 (S11E38 L=278)   325    119        0 ( 421/ 761- 514/ 844) -51.3 E41-E35 S13-S08
      6  430  620 (S22E45 L=271)   585    119        2 ( 401/ 556- 461/ 676)  41.2 E49-E41 S25-S18
      7 1585  623 (S30W38 L=354)   349     84        0 (1550/ 597-1619/ 656)  31.0 W35-W40 S31-S28
      8  775  529 (S31E21 L=295)   225     58        1 ( 751/ 506- 796/ 552)  36.9 E23-E19 S32-S29
      9 1731  520 (S38W59 L= 15)   440     75        2 (1689/ 506-1767/ 540)  -4.8 W53-W65 S39-S36
     10  746  374 (S41E29 L=287)   402     52        1 ( 725/ 344- 769/ 402)  63.4 E30-E26 S43-S39
     11  677  275 (S49E43 L=274)   528     85        1 ( 645/ 244- 717/ 299)  17.1 E49-E36 S51-S46
  

    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