Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  696 1359 (N20E19 L=106)   278     83        0 ( 663/1311- 736/1402)  50.2 E21-E16 N16-N22
      2  524 1287 (N16E31 L= 94)   387    110        0 ( 473/1256- 562/1309) -26.6 E34-E28 N14-N17
      3  206 1254 (N17E59 L= 67)  1034    282        3 ( 146/1199- 259/1303) -29.7 E66-E52 N12-N20
      4 1217 1162 (N05W12 L=137)   237     61        0 (1182/1129-1245/1188) -45.0 W10-W13 N02-N06
      5 1520 1019 (S05W31 L=156)   257     62        0 (1500/ 984-1538/1054) -53.1 W29-W32 S06-S02
      6 1696  951 (S09W44 L=169)   782    204        0 (1660/ 844-1718/1053) -72.9 W41-W46 S15-S02
      7  123  947 (S01E70 L= 55)   850    224        0 (  80/ 922- 167/ 971)  -7.1 E80-E63 S02-N00
      8  385  712 (S17E45 L= 80)   368     97        2 ( 357/ 686- 415/ 745) -26.6 E47-E43 S18-S15
      9  502  656 (S22E37 L= 88)   491     92        1 ( 478/ 607- 523/ 703) -90.0 E38-E35 S24-S18
     10 1845  665 (S27W62 L=187)   545     35        0 (1828/ 642-1862/ 688)  45.0 W60-W63 S28-S25
     11  388  605 (S24E48 L= 77)   327    101        0 ( 362/ 580- 431/ 650) -45.0 E49-E44 S25-S21
     12 1103  492 (S36W03 L=128)   677    188        0 (1020/ 470-1204/ 511)  -4.4 E03-W10 S37-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