Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1441 1473 (N18W39 L=327)   242     83        1 (1408/1445-1469/1500)   9.5 W35-W41 N17-N19
      2  611 1388 (N37E17 L=272)   241     64        0 ( 580/1373- 647/1402)  -9.5 E19-E13 N36-N38
      3 1226 1348 (N18W21 L=309)   197     76        0 (1200/1332-1267/1375) -63.4 W19-W22 N15-N19
      4 1251 1153 (N06W16 L=304)   906    264        3 (1123/1082-1375/1182) -25.0 W08-W23 N00-N10
      5  825 1070 (N13E10 L=279)   393    167        3 ( 739/1021- 896/1114)   0.0 E16-E05 N11-N14
      6 1797  950 (S21W48 L=336)   348     74        1 (1781/ 906-1815/ 983)  31.0 W45-W50 S23-S19
      7  192  880 (N16E62 L=227)  2473    476        4 (  90/ 835- 307/ 929) -16.9 E74-E50 N11-N20
      8 1366  663 (S24W10 L=299)  1474    344        4 (1312/ 514-1430/ 801)  49.1 W03-W17 S33-S17
      9  468  644 (S02E44 L=244)   256     89        1 ( 425/ 617- 496/ 665) -36.9 E47-E42 S04-N00
     10  756  516 (S16E30 L=258)   262     41        0 ( 730/ 501- 783/ 534) -33.7 E31-E28 S17-S15
     11  639  418 (S20E42 L=246)   466    157        2 ( 580/ 401- 724/ 434) -24.0 E46-E37 S22-S17
  

    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