Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  767 1877 (N41E50 L=178)   975    250        0 ( 750/1829- 786/1931) -35.5 E58-E44 N36-N47
      2  594 1507 (N11E39 L=190)  1214    338        0 ( 532/1367- 676/1649) -87.1 E41-E36 N01-N21
      3 1368 1216 (N12W15 L=243)  2030    487        1 (1275/1009-1519/1440) -38.4 W04-W28 N04-N23
      4 1521 1403 (N28W21 L=249)   295     38        0 (1494/1377-1543/1432) -26.6 W18-W22 N26-N29
      5  717 1315 (N02E25 L=204)  1475    398        2 ( 530/1219- 848/1358)   2.6 E37-E15 S01-N05
      6  277 1284 (S08E56 L=172)   555    111        1 ( 241/1237- 322/1330)  66.8 E58-E52 S11-S04
      7  839 1122 (S06E13 L=216)   591    161        0 ( 813/1044- 858/1205) -60.9 E15-E10 S10-S01
      8 1821 1131 (N22W49 L=278)   893    122        2 (1793/1065-1838/1189) -50.2 W45-W51 N17-N25
      9  669 1068 (S13E22 L=207)  1097    254        5 ( 569/1021- 779/1121)  39.3 E27-E15 S17-S08
     10 1822  876 (N07W56 L=285)  1910    290        1 (1788/ 757-1864/ 993) -34.5 W49-W66 N02-N13
     11 1054  700 (S24W10 L=239)  1121    386        1 ( 867/ 671-1245/ 739)  27.6 W00-W22 S30-S18
     12  495  699 (S39E28 L=201)  1649    422        5 ( 345/ 679- 667/ 725)  16.5 E41-E14 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