Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  476 1729 (N41E56 L=295)   481     46        1 ( 447/1706- 502/1743)  -7.1 E61-E52 N39-N42
      2 1733 1404 (N24W51 L= 42)   367     50        0 (1712/1374-1756/1431) -56.3 W48-W52 N22-N25
      3  206  817 (S18E61 L=289)   528    105        0 ( 192/ 766- 224/ 865)  50.2 E64-E58 S21-S15
      4  348  755 (S22E47 L=304)   378    121        1 ( 304/ 722- 380/ 809)  36.9 E52-E44 S24-S18
      5  455  766 (S22E38 L=313)   423     67        2 ( 431/ 732- 479/ 801) -76.0 E39-E36 S23-S19
      6  399  514 (S38E51 L=300)  1039     93        4 ( 362/ 482- 445/ 549)  17.1 E57-E44 S40-S35
      7 1374  491 (S36W30 L= 21)   328     78        0 (1360/ 449-1393/ 538) -56.3 W27-W32 S38-S32
      8 1487  512 (S33W38 L= 29)   255     36        1 (1464/ 495-1507/ 532)  18.4 W36-W40 S34-S31
  

    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