Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  653 1772 (N33E51 L=261)   458    139        0 ( 630/1735- 671/1800) -26.6 E55-E47 N30-N35
      2  385 1523 (N11E58 L=254)  5686    831        9 ( 297/1308- 496/1748) -71.6 E66-E53 S03-N27
      3  461 1498 (N11E50 L=263)   603    130        2 ( 429/1456- 492/1538) -71.6 E53-E47 N08-N14
      4  128 1217 (S11E76 L=237)  1248    120        2 ( 112/1179- 142/1253) -53.1 E79-E72 S13-S09
      5  431 1088 (S13E39 L=274)   612    101        1 ( 402/1031- 468/1130)  81.9 E40-E35 S17-S10
      6  416  824 (S30E36 L=277)   327     76        0 ( 397/ 778- 437/ 855)  68.2 E36-E33 S32-S27
      7  971  765 (S19W03 L=316)   521    181        0 ( 896/ 741-1081/ 781)  33.7 W00-W09 S22-S16
      8  294  684 (S42E47 L=265)  1131    299        1 ( 228/ 638- 366/ 728)  29.4 E55-E39 S46-S37
      9 1672  626 (S09W54 L=  6)   427     66        0 (1643/ 596-1700/ 659)  68.2 W52-W55 S11-S06
     10 1712  574 (S11W61 L= 13)   884    205        0 (1656/ 548-1765/ 594)   8.7 W54-W67 S12-S08
     11 1387  449 (S27W41 L=353)   635    206        0 (1299/ 425-1462/ 479)  13.0 W34-W47 S28-S24
     12 1351  268 (S39W51 L=  3)  1066    246        1 (1328/ 215-1376/ 341) -32.5 W45-W57 S42-S34
  

    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