Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  243 1333 (N02E64 L= 49)   586    102        0 ( 206/1310- 272/1362)  26.6 E68-E60 N00-N03
      2  959 1295 (N13E09 L=103)   363    127        0 ( 921/1258-1008/1360) -26.6 E13-E05 N11-N16
      3  376 1044 (S13E43 L= 70)   479    139        0 ( 358/ 965- 393/1107) -77.5 E44-E41 S18-S09
      4 1117  954 (S04W07 L=119)   335     85        1 (1093/ 908-1135/1004) -80.5 W05-W08 S06-N00
      5  933  937 (S09E04 L=109)   370    118        2 ( 901/ 868- 959/ 997) -82.9 E05-E02 S13-S05
      6  791  792 (S20E10 L=103)   476    207        1 ( 693/ 754- 883/ 841)  39.3 E15-E04 S24-S15
      7 1075  738 (S17W09 L=122)   220    100        1 (1041/ 704-1121/ 757)  45.0 W07-W11 S20-S16
      8 1144  684 (S19W15 L=127)   298    134        1 (1101/ 631-1197/ 727)  76.0 W12-W17 S23-S15
      9  770  671 (S28E09 L=104)   344    156        0 ( 701/ 624- 829/ 723)  52.1 E12-E05 S32-S23
     10 1269  689 (S16W22 L=135)   215     67        0 (1247/ 663-1307/ 712)  63.4 W20-W23 S18-S14
     11 1345  322 (S38W44 L=156)   930    181        1 (1266/ 302-1421/ 345)  14.0 W36-W48 S39-S34
     12 1395  255 (S41W53 L=166)   455     64        1 (1369/ 237-1417/ 271)  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