Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  961 1298 (N13E09 L=103)   480    134        0 ( 920/1257-1016/1366) -26.6 E13-E05 N11-N16
      2  245 1332 (N02E63 L= 49)   586    110        1 ( 209/1310- 272/1365)  36.9 E67-E60 N00-N03
      3  827 1121 (S00E13 L= 99)   290     50        0 ( 808/1088- 845/1152) -26.6 E15-E11 S01-N01
      4  379 1040 (S13E42 L= 70)   579    134        0 ( 361/ 964- 397/1104) -77.5 E44-E41 S18-S09
      5 1121  958 (S04W07 L=119)   376     96        1 (1094/ 907-1140/1009) -80.5 W05-W08 S06-N00
      6  940  936 (S08E03 L=109)   435    125        3 ( 913/ 874- 966/ 993) -81.9 E04-E01 S12-S05
      7  791  789 (S20E09 L=103)   628    223        0 ( 679/ 753- 890/ 840)  36.9 E16-E04 S24-S15
      8  781  674 (S28E08 L=104)   515    199        0 ( 706/ 622- 871/ 730)  48.0 E12-E03 S32-S22
      9 1151  683 (S19W15 L=127)   418    144        0 (1106/ 627-1210/ 726)  66.8 W13-W18 S23-S16
     10 1274  689 (S16W22 L=135)   250     59        1 (1250/ 665-1311/ 710)  45.0 W21-W24 S17-S14
     11 1347  321 (S38W44 L=156)  1007    181        1 (1269/ 303-1422/ 346)  24.4 W37-W48 S40-S34
     12 1394  253 (S42W53 L=165)   524     68        1 (1367/ 234-1420/ 270)  90.0 W50-W55 S43-S39
  

    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