Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1069 1505 (N32E03 L=287)   319     56        1 (1054/1467-1086/1543)  90.0 E03-E01 N29-N35
      2 1324 1190 (N14W17 L=307)   290    111        0 (1300/1138-1342/1238) -63.4 W15-W18 N11-N17
      3 1646 1179 (N17W41 L=331)   317     75        1 (1612/1140-1675/1205) -31.0 W38-W43 N14-N18
      4 1874 1063 (N11W65 L=355)  1758    217        2 (1853/ 959-1888/1171)  85.6 W62-W68 N05-N18
      5  144 1008 (S10E70 L=220)   832    118        0 ( 101/ 977- 196/1037)  -5.4 E82-E61 S11-S08
      6  823  884 (S09E11 L=279)   636    180        1 ( 771/ 809- 864/ 958)  68.2 E13-E08 S14-S04
      7  323  878 (S15E48 L=242)   804    108        2 ( 266/ 843- 368/ 924)   5.7 E53-E43 S17-S12
      8  210  794 (S22E63 L=227)  2196    192        6 ( 166/ 716- 262/ 881) -66.0 E70-E56 S27-S16
      9  300  748 (S24E52 L=238)   867    105        2 ( 273/ 699- 320/ 802) -81.9 E54-E49 S27-S20
     10 1549  575 (S22W43 L=333)   510    166        0 (1484/ 545-1597/ 625)  26.6 W37-W45 S23-S17
     11 1268  406 (S37W27 L=317)   410     98        0 (1257/ 357-1278/ 453) -63.4 W24-W29 S40-S32
  

    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