Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1104 1848 (N55E01 L=284)   345     49        1 (1081/1822-1126/1870)  63.4 E03-W00 N51-N57
      2  459 1625 (N31E50 L=236)   444    103        2 ( 422/1589- 497/1664)  45.0 E52-E47 N27-N33
      3  925 1563 (N28E11 L=275)   377    154        0 ( 860/1530- 992/1577)  20.6 E14-E06 N25-N29
      4 1525 1448 (N26W31 L=317)   239     78        0 (1489/1426-1565/1469) -18.4 W28-W34 N24-N26
      5  825 1134 (N00E13 L=273)   205     48        0 ( 802/1105- 853/1158)  45.0 E13-E10 S01-N02
      6  756  706 (S26E15 L=270)   238     74        0 ( 724/ 677- 796/ 727)  31.0 E17-E12 S27-S24
      7  886  497 (S39E06 L=280)   307     67        0 ( 852/ 483- 928/ 513)   0.0 E08-E02 S39-S37
      8 1135  462 (S39W13 L=299)  1064    178        3 (1058/ 421-1191/ 508)  11.3 W07-W17 S42-S35
      9 1594  355 (S41W62 L=348)   869    115        0 (1566/ 331-1650/ 399)  -6.7 W55-W72 S42-S38
  

    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